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A METHOD FOR COMPUTING LUMINOUS EFFICIENCIES 
FROM METEOR DATA 

By Barrett S. Baldwin, Jr-, and H- Julian Allen 
Ames Research Center 


SUMMARY 


The details of a model for the ablation of a stone meteoroid passing 
through the atmosphere are presented. The model allows for a decrease In 
mean density due to frothing, which results from the intense heating at low 
atmospheric pressures- The model is applied to an analysis of meteor data to 
evaluate the luminous efficiency factor- The present results are a refine- 
ment of a previously published analysis since new information on the proper- 
ties of meteoritic materials has been utilized. A method has also been found 
for taking into account a delay in light production resulting from the time 
required to vaporize chunks of froth in the meteor wake- 

The introduction of a frothing ablation model in an analysis of super- 
Schmidt meteor data allows an interpretation in terms of originally solid 
meteoroids for meteors previously thought to have been produced by fragile 
low-density meteoroids. At the same time, the computed luminous efficiency 
factor, t 0 p y is found to vary substantially with velocity. Above 35 km/sec 
the value of T 0 p is close to that recommended by Verniani . As the velocity 
is reduced below 25 km/sec r G p increases to a level about six times greater. 


INTRODUCTION 


The fragmentation concept was originated in 1955 by Jacchia (ref- l) as 
an explanation of the deceleration anomaly of faint meteors. The idea was 
supported by Whipple* s earlier conclusion that cometary meteoroids were prob- 
ably porous and fragile (see ref- 2 ). This view was widely accepted and 
attempts have been made to predict the behavior of such objects in the 
atmosphere, for example, Opik*s dust ball theory (ref. 3 ) • 

An alternative explanation of the observed anomalous decelerations and 
low densities was presented in reference 4. The explanation is based on 
results of ground-based arc- jet tests where it was shown that many materials, 
including those found in meteorites, develop low-density froths under condi- 
tions of high heating rates and low pressures- Frothing can result either 
from boiling of the melted material or expansion of volatile components con- 
tained in the melt. It can be concluded from these tests that a sample from 
an ordinary stone meteorite weighing one gram or less, if impelled into the 
upper atmosphere at meteoric velocity, would develop a shell of low-density 
froth. The mean density would thereby be lowered and could decrease during 
the trajectory, leading to an anomalously large deceleration. Although this 


finding does not preclude the possibility that low-density porous meteoroids 
exist in space , it weakens some of the arguments that originally led to that 
hypothesis . 

The view that a majority of meteoroids may be solid rather than porous, 
suggested in reference 4, was based in part on a mathematical model for the 
behavior of such meteoroids during ablation. The model was applied to an 
analysis of super-Schmidt meteor data to determine values of the luminous 
efficiency factor r Q p and the mean meteoroid density p m . Previous meteor 
analyses that do not include a model for the ablation process can only predict 
the value of the ratio T op /pj|* The level of T op computed from data con- 
tained in reference 5 compared favorably with previous estimates. However, in 
contrast to the prevailing assumption (see, e-g., refs. 5 to 8) the luminous 
efficiency factor was found to increase with decreasing velocity rather than 
remaining constant over the entire range of meteor velocities. The predicted 
velocity dependence of the meteoroid flux in space is modified if T op is 
variable rather than constant. It was shown in reference 4 that the resulting 
modification leads to a better agreement with velocities predicted from 
overlapping radar data. 

The present report explains the details of the mathematical model used 
for analyzing meteor data in reference 4. It also improves the model by 
including laboratory measurements of the properties of meteoritic materials. 

A method has also been found for taking into account a delay in light produc- 
tion when froth is sloughed off in relatively large pieces. It is shown that 
this delay is partially responsible for apparent increases of luminous effi- 
ciency along a trajectory. The evidence in support of the frothing- sloughing 
ablation model is thereby strengthened. 


SYMBOLS 


A meteoroid frontal area 

A p froth particle frontal area 

A rn geometric parameter in ablation model (eq. (A10) ) 

A rs geometric parameter in ablation model (eq. (A2) ) 

A 0 frontal area of solid inner core 

O 

B fraction of mass ablated as froth (eq. ( 19 )) 

C heat -transfer parameter (eq. (C 3 )) 

C D drag coefficient 

C p skin-friction coefficient 
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C H 

^Hcvu 

C Heu 

C Huu 

CHu 

c lJ C 2 

cos Zp 

F 

G 

glj>g2;g3 

h 

i 

k 

l f 

l s 

M 

m 

m 

m i 

m l 


uncorrected skin-friction coefficient 

aerodynamic heat -transfer coefficient 

convective heat-transfer coefficient 

equilibrium radiative heat-transfer coefficient 

nonequilibrium radiative heat-transfer coefficient 

uncorrected heat-transfer coefficient 

parameter in viscosity law (eq. (Bll)) 

parameters in vapor pressure law (eq. (Oil)) 

cosine of angle between meteor path and the vertical 

factor (eq. (D8)) 

correction factor to account for nonlinear velocity profile 

factors to allow for nonlinearities and nonuniformities in the 
viscous flow of the melt 

instantaneous meteor altitide 

photographic luminous intensity 

parameter in viscosity law (eq. (Bll)) 

thermal conductivity 

length of cylindrical froth shell (fig- l8) 
length of cylindrical solid core (fig- l8) 
molecular weight of vapor 
meteoroid mass 

rate of ablation of meteoroid mass with time 

meteoroid mass based on T 0 pq 

mass of froth particle in wake 

mass of liquid drop in wake 

mass of solid inner core 



m s 

P v 

P v ^Py( T ) 



a 


AL 


€ 


s 


rate of melting of solid core with time 
heat -transfer parameter (eq. ( C3 ) ) 
equilibrium vapor pressure 
meteoroid total volume 
volume of solid inner core 
ratio (eq. (l2)) 

radius of meteoroid frontal area 
radius of froth particle in wake 
radius of liquid drop in wake 
nose radius 

radius of frontal area of solid core 
shape factor (eq. (4)) 
shape factor of hemisphere^ 
meteoroid surface temperature at stagnation point 
temperature of fusion 
time 

tangential velocity of melt 

tangential surface velocity of melt at control surface 
meteoroid velocity 
meteoroid acceleration 

rate of recession of solid surface at stagnation point 

heat-transfer parameter (eq. (C3)) 

lag distance of froth in meteoroid wake (eq. (22)) 

thickness of melted layer at control surface 

surface emissivity 

specific heat of ablation 



k 



specific heat of fusion (latent plus stored) 
specific heat of vaporization (latent plus stored) 

9 intensity lag time 

time to ahlate froth particle to nothing in wake 
9 ^ time to vaporize liquid drop 

jjl viscosity of solid 

\±£ viscosity of froth 

p air density 

*p ratio of air density to sea-level air density 

p Q sea-level air density 

p 21 density ratio across normal shock 

p m meteoroid mean density 

p m rate of change with time of meteoroid mean density 
p mf froth density 

P ms solid density 

a Stefan-Boltzmann constant 

a 1 heat-transfer parameter (eq. (C3)) 

T 0 p luminous efficiency factor (zero mag-cgs) 

T o p i assigned constant luminous efficiency factor (6.46x10” 19 zero mag-cgs) 
T opc corrected luminous efficiency factor 
t luminous efficiency (dimensionless) 

ir 

T g surface tension 

to rate of vaporization per unit area 

to total rate of ablation per unit area due to evaporation and explosive 

spalling 

to^ rate of vaporization per unit area in a viscous boundary layer 
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rate of vaporization per unit area in a vacuum 

a) rate of melting of solid core mass per unit area and time (eq. (Al4)) 
s 

Sub scripts 

f final point on trajectory (also froth) 
i integer denoting a point on a meteor trajectory 

j integer representing the number of iterations 

-oo prior to entry into the atmosphere 


AKA.LYSIS 


In this section a set of relations will be exhibited that can be used to 
derive values of the luminous efficiency factor T op from meteor data. These 
relations are, in part, based on an ablation model which takes into account 
the production and viscous flow of froth in the heated region at the front of 
a meteoroid. 


Basic Meteor Equations 

In reference 5, meteor data, reduced from photographic data by the 
methods described in references 9 and 10, are presented as values at several 
points on the trajectories of each meteor of the following quantities: 

h altitude 

m 1 mass 

V velocity 

V acceleration 

An identifying trail number and a value of cos Z^, the cosine of the angle 
between the meteor trajectory and the vertical, are given for each meteor. 
Much additional information, not of direct concern here, is included. The 
values of mass are based on an assumed constant value of the luminous 

efficiency factor 
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T op i = 6.46x10 


~ 19 (zero mag-cgs ) 1 

and a final mass equal to zero. 

The luminosity and dynamic equations used in reference 5 are 

Ip = - 75 T 0 pV 3 ~ (luminosity equation ) 1 

and 

v = - ~ Cr> 

2 v m 


V = - % C D | pv 2 


(dynamic equation) 


( 1 ) 

( 2 ) 


At the outset, we shall take the masses in these two relations to be equal as 
in references 2 to 8 . However, an alternative assumption allowing for a lag 
in the time required for light to be produced from froth particles in the wake 
is investigated and found to have merit. 


The values of mass m ± in reference 5 are based on the assumed luminous 
efficiency factor T 0 pq and in accord with an integral of equation (i) . The 
values of mass can be converted to correspond to a different luminous effi- 
ciency factor T 0 p by integration backward along the trajectory using the 
relation 


m 


r m i T opi 

Jo T op 


dm x 


(3) 


derived from equation (l). If T 0 p is assumed to be constant between points 
on the trajectory, the integral in equation (3) can be replaced by a sum. 

As in references 4 and 7 we shall assume that the shape of the meteoroid 
remains similar as it ablates during entry. Then the shape factor will be a 
constant and is expressed as 


s F = a/q 2/3 (4) 

where Q is the meteoroid volume given by 

Q = m/p m (5) 


^The units used in this report for all equations except the luminosity 
equation (l) are MKS. In the discussions, other units, such as grams and 
kilometers per second, may be specified. The units in the luminosity equa- 
tion are those of reference 5 - The photographic intensity Ip is expressed 
in units of the intensity of a zero magnitude star. The unit of the luminous 
efficiency factor T 0 p is zero magnitude sec 4 gram“ 1 cm“ 3 , abbreviated as 
zero mag-cgs. For use in other equations, the mass determined from equa- 
tion (l) is multiplied by 1CT 3 to convert to MKS units. Methods for express- 
ing the luminous efficiency in other forms that appear in the literature are 
considered in a later section. 
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in terms of the mass m and the meteoroid mean density p m . From equations 
(2), (4), and (5) the mean density of the meteoroid can he expressed as 


j/ s F c D pv 2 \ 3 /2 

P m •Jm\ 2 V ) 


(6) 


In this report, values of the air density p corresponding to the altitudes 
given in the data of reference 5 'will he taken from reference 11. The meteor- 
oid front face will he assumed to he a segment of a sphere and the drag 
coefficient Cp evaluated according to reference 12. 

It can he seen that application of equations (3) and (6) to the reduced 
meteor data of reference 5 leads to evaluation of^the ratio SpT 0 p/Pm each 
point on the trajectories where values of h, V, V, and m x are given. Fur- 
thermore, if values for the shape factor Sp and meteoroid density Pm are 
assumed as in reference 6, the luminous efficiency factor T 0 p can he evalu- 
ated. Alternatively, an additional relation involving these quantities can he 
arrived at hy consideration of the ablation process occurring during atmo- 
spheric flight as in references 7 and 13. The ablation model developed for 
this purpose in the present report differs from those previously employed in 
that it takes into account the frothing and sloughing of meteoritic materials 
as observed in arc- jet tests. Application of this model to reduced meteor 
data from reference 5 leads to evaluation of t op an( i Pm individually as well 
as the ratio T op /p^ resulting from previous analyses. For that purpose it 
remains necessary to assume a value for the shape factor Sp, hut the effect 
on the results of variations in this quantity is investigated. 



Ablation Model 

The arc -jet tests on which our 
model is based are described in ref- 
erences 4 and l4. Figure 1 illus- 
trates the processes observed in the 
tests, which we assume also occur in 
meteor flight. The solid stone core 
melts at the front face and runs to 
the side with only a small amount of 
vaporization. Upon reaching a 
region of low enough pressure, the 
melted stone froths. The frothing 
can he caused hy volatile components 
in the melt, or, in the absence of 
such components, hy boiling. After 

reaching the sides, where the heating from the air is less severe, the froth 
solidifies. Chunks of froth intermittently break off and, in the case of a 
meteor, are assumed to undergo further ablation and eventually vaporize in the 
wake. In a part of the range of test conditions covered, the froth ran off 
the model in the liquid form rather than solidifying. This occurred at condi- 
tions corresponding to the later parts of slow meteor trajectories. 


Figure 1. - Sketch of the cylinder (solid) 

hemisphere (froth) assumed in the meteoroid 
ablation model. 


8 



Since the tests were run under conditions that did not correspond com- 
pletely to meteor flight, a discussion of the differences is in order. The 

heating loads and impact pressures of the tests corresponded to slow 

(l5 km/sec) meteors at altitudes from 70 to 100 km. The effective altitude 
was held constant in each test, and acceleration effects were not simulated. 
Models made of material taken from stone and iron meteorites were tested and 
found to react in the same manner as similar terrestrial materials. The 
models were originally round-nosed cylinders with a radius of 1 cm, weighing 
about 30 grams. Tests as long as 115 seconds were required at the higher 
effective altitudes to achieve the conversion of an appreciable fraction of 
the material to froth. It is important to note that frothing is a surface 
phenomenon rather than one that occurs throughout the volume of the solid mate- 
rial. Therefore, frothing of the same fraction of mass will occur in a 

shorter time in a smaller object. For example, if the time required for melt- 
ing a given fraction of the mass is taken to be proportional to the ratio of 
volume to surface area, the same degree of melting would occur in a 1 gram 
mass in one-third the time required for a 27“gram mass. Therefore, the frac- 
tion of a 30 -gram model converted to froth in 115 seconds at an altitude of 
98 km would correspond to a test time of about 36 seconds for a 1-gram mass at 
the same altitude. Since the heat load is approximately proportional to air 
density, the same degree of frothing would be expected in 3-6 seconds in a 
1-gram mass at an altitude of 86 km where the air density is about 10 times 
greater than at 98 km. This forms a part of the basis of our assumption that 
meteoroids in the mass range below 1 gram will undergo a degree of frothing 
comparable to that observed in the tests. 

The effective velocity of the tests (about 15 km/ sec) corresponds to the 
velocity of slow meteors. Since the heating is approximately proportional to 
the cube of velocity, an increase in velocity by a factor of 3 would decrease 
the time required for melting at a given altitude by about a factor of 27* 
Therefore, a fast meteor below 1 gram in mass receives sufficient heat to 
froth at altitudes above 98 km. It is possible that other effects not 
observed in the tests would occur at velocities above those reached. For 
example, vaporization might take place rapidly enough to prevent the degree of 
froth build-up observed in the tests. We do not have direct evidence to rule 
out this possibility. However, the calculations to be discussed later in this 
report indicate that for all of the meteors analyzed most of the solid mate- 
rial at the front face is removed by viscous flow of the melted material, 
rather than by vaporization. Our conclusion from this is that the eventual 
vaporization of meteoroids of this size takes place largely in the wake rather 
than on the parent meteoroid. An important link in this reasoning is the 
finding that, at the temperatures involved, the viscosity of meteoritic stone 
(and similar terrestrial materials) is several orders of magnitude lower than 
that of a glass with high silicate content. 

The major assumptions utilized in the ablation model are as follows: 

1. The meteoroid front face, comprised partly of froth and partly of a 
thin layer of melted material covering the solid core, remains shaped as the 
surface of a segment of a sphere with the radius of the frontal area r a 
constant fraction of the nose radius r n during ablation. 
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2 . The ratio of solid volume to total volume Q s /Q is proportional to 
the cube of the ratio of the frontal area radii r s /r (see fig. l) . 

3. During ablation the spherical front face of the solid remains cen- 
tered in the froth and the entire front face remains oriented -with respect to 
the air stream as shown in figure 1 . 

4 . Ablation of the solid material takes place by vaporization and vis- 
cous flow with negligible effects from surface undulations such as ring waves 
(ref. 15). The solid core does not fragment. 

5 . The solid core is composed of a meteoritic material for which the 
properties such as thermal heat conductivity , viscosity, etc., are known. 

The origin and implications of the assumptions are discussed in appen- 
dix A. The viscous flow of melted material resulting from the impinging air- 
stream is analyzed only for the thin layer covering the front of the solid 
core. The ablation model based on the foregoing assumptions is thus in a 
sense incomplete in that no attempt is made to predict the degree of froth 
build-up or to follow the mode of ablation of the froth (whether by vaporiza- 
tion, viscous flow, or fragmentation). Instead, the degree of froth accumula- 
tion is deduced from the meteor data with the aid of the foregoing assumptions 
plus those implicit in the basic meteor equations. The geometric assumptions 
are made in lieu of an attempt to guess the original shape and predict 
evolutions of shape from a detailed study of the ablation processes occurring. 


An expression for the rate of melting of solid core mass per unit area 
and time is derived in appendix A in the form 


u s = 


Af s 


Sp 


Pm 

m 


2/3 


Pms 

Pm 


Pm - PmfV /3 


ms 


“ Pmf 


m + 


mp ms p mf 


( p m p mf) ( p ms 




p m 


m 


( 7 ) 


where 

A rs geometric parameter that determines the radius of the solid frontal area 
relative to the radius of the total frontal area 

Sp shape factor 

m meteoroid mass 

p m mean meteoroid density 

p mf dens i' t y ° f froth 

p ms ( 3- ens i' t y of solid core 

The rate of melting of solid mass per unit area is related to the values 
and derivatives of the meteoroid mass and mean density by this equation when 
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Control surface 




Figure 2 . - Coordinate system for the viscous 
flow of the melt in the thin layer at the 
front of the solid core. 


the parameters A rs , Sp, p^g, and 
Pmf are assigned. The values of 
the parameters used for most of the 
calculations are listed in table I, 
but the effects on the results of 
reassignment of the values will be 
considered . 

In appendix B, equations are 
derived that express the conserva- 
tion of mass, momentum, and energy 
for the viscous flow of the thin 
layer of melted material at the 
front of the solid core . Figure 2 
indicates that part of the meteoroid 
under consideration. The relations 
derived in appendix B are combined 
to form the equation 


_k \ 2 p ms 
k) M-(T) 


(T - T f ) 2 G(T) 



s/TtCppV^ 


( 8 ) 


where 

k 

£f 

p ms 

m-(t) 

T 


T f 

G(T) 



m 


^m 


oo 


thermal conductivity of melted meteoroid 
specific heat of fusion plus sensible heat 
density of solid meteoroid core 
viscosity of melted meteoroid 
stagnation-point surface temperature 
temperature of fusion of meteoroid 

correction factor for nonlinear velocity profile (appendix B) 

ratio of nose radius to radius of frontal area 

shape factor 

meteoroid mass 

meteoroid mean density 

rate of melting of the meteoroid solid core per unit area and time 

(eq. (7)) 

rate of evaporation of solid core material per unit area and time 
(eq. (C12)) 
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Cp tangential surface force coefficient 
p air density 
V meteor velocity 

The material property laws assumed and the aerodynamic laws utilized are dis- 
cussed in appendixes B and C. Most of the results in this report are Based on 
material properties appropriate for carbonaceous chondrites- However, the 
effect on the results of variations in properties is investigated. 

If the quantities appearing in the auxiliary relations are considered 
evaluated, substituting equation (7) into (8) would lead to a complicated 
expression which can be indicated symbolically as 

h,V,V) = 0 (9) 

This ablation equation provides an additional relation between the quantities 
appearing in the basic meteor equations. It can be seen that insertion of 
values of h, V, and V from reduced meteor data into equation (6), a form of 
the dynamic equation, and the ablation equation (9) leads to evaluation of m 
and p m along the meteor trajectory. That information plus the values of m^ 
from the reduced meteor data can be used in equation ( 3 ) 9 a form of the 
luminosity equation, to determine the luminous efficiency factor Top • As a 
result of this procedure, values of both T 0 p and p m are obtained (only the 
ratio T 0 p/p^ is found from the use of the luminosity and dynamic equations 
without an ablation equation) . 


Method of Solution 


An iterative procedure is used for the calculation of values of Top and 
Pm ea ch point on the trajectory of a meteor where values of h, V, V, and 
in-L are given in reference 5* A computer program has been developed for that 
purpose. The computations performed in each stage of the iteration are 
indicated in the following numbered paragraphs . 


1. The mass m along the trajectory is assigned the values from the 
previous iteration, and equation (6) is used to compute at each 

trajectory point. 


The evaluation of m is based on the relation 


„ dm _ „ d ln(m) 

m = -V cos Z R — = -V cos Z R m dh 


obtained with the aid of the trajectory relation 

dh 
dt 


dh ,r n 

— = -V COS Z R 


( 10 ) 


The logarithmic differentiation is employed because of the relative constancy 
of dZn(m)/dh compared to dm/dh along the trajectory. Similarly, the 
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evaluation of is based on 


Pm = -V cos Z rP] 


d InCpm) 


m 


dh 


( 11 ) 


The derivative dZn(m)/dh is approximated at midpoints between trajec- 
tory points by simple differences according to the relation 


d Xn(m) 


dh 




ln(m i+1 ) - Zn(m i ) 
h i+i - h i 


where the subscript i denotes the ith trajectory point. A linear inter- 
polation is used to compute derivatives at the trajectory points. 


d Zn(m) 

f h i+i “ h i \ 

[~d Zn(m) 

. ~ h i-i 'N 

d 7n(m) 


dh 

j_ vH+i " v 

| dh 

_ l v'i+i " h i -x) 

± p 

; dh 

1 

1+ 2 


The derivative dXn(p m )/dh is evaluated by the same method. This method is 
not applicable at the first and last trajectory points, which are therefore 
dropped from the subsequent calculations. For that reason, only meteors with 
data given for at least three trajectory points can be used in our machine 
program in its present form. 

2. The quantities m and along the trajectory are evaluated 
(eqs. (lO) and (ll)) using the numerical procedure described above. 

3* Values of w s are computed by means of equation (7) • 

4. The stagnation-point temperature T and the quantities p(T), g(t), 
£3, and Cp in equation (8) are computed with the aid of auxiliary relations 
given in appendixes A, B, and C. 


The viscous flow equation (8) is redundant in terms of the quantities 
computed in the previous steps. Therefore , it can be used to readjust the 
values of mass that were assigned at the beginning. For that purpose the 
ratio of the two sides of equation (8) 


R = 


a o 1 / 2 /"mV-/ 3 p / ~\ 

A rn S F " w ) 

^C F pV 2 (j~y ( T _ Tf ) 2 G(T) 


(12) 


serves as an indicator of the needed change. By trial and error it has been 
found that if R is not equal to 1, as it should be, readjusting the luminous 
efficiency factor according to the relation 


T °Pj 


= R l/2 ^ 


r °Pj-i 


(13) 
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will give a value of R closer to 1 in the next iteration. The subscripts 
j and j-1 refer to the jth and (j-l)th iterations, respectively. 


5* The next step in the calculation is the evaluation of R and r Q p at 
each point of the trajectory (except the first and last points) using 
equations (12) and ( 13 ) - 


6. Equation ( 3 ) can now be used in the form 

T °pj-i 


m 




1 op . 




( 1*0 


to compute the adjusted values of mass along the trajectory. The subscript j 
refers to the stage of the iteration and the subscript i to the point on the 
trajectory. For the final point on the trajectory, the relation 


m 


jf 


' T ^A 

r °?j It- 1 


j- 1 ,? 


(15) 


is used, where the subscript f-1 refers to the next to last point. The 
luminous efficiency factor for the last point is not computed because of 
failure at that point of the method of numerical differentiation employed in 
step 2 . 

After completion of step 6, new values of mass along the trajectory are 
available for use in the next iteration starting with step 1. For most of the 
reduced meteor data from reference 5 treated in this report , 10 iterations 
lead to values of R that differ from 1 by less than 10 percent. It has been 
found that the resulting values of m, p^, and T 0 p along the trajectory are 
independent of the luminous efficiency factor assumed at the beginning of the 
iteration to provide starting values of the mass. 

In the discussion following equation ( 9 ) , it was stated that values of 
mass m and mean meteoroid density could be derived from the meteor data 

with the aid of the dynamic and ablation equations without reference to the 
luminosity equation. That purpose can be accomplished in the foregoing itera- 
tion procedure by merely declining to identify the quantity r Q p computed in 
equation ( 13 ) as the luminous efficiency factor. In that case, equation (l4) 
can be regarded as one of the mathematical steps in the iteration procedure 
rather than identifying it with the luminosity equation. Therefore, if modi- 
fications of the luminosity equation are considered, it is not necessary to 
change equation (l4) or any part of the iteration procedure. The values of 
mass and meteoroid density computed from the iteration are thereby determined 
from the dynamic and ablation equations alone. They can therefore be used in 
a modified form of the luminosity equation to determine a luminous efficiency 
factor T 0 p C that is not equal to the dummy variable T 0 p used in the itera- 
tion. We shall have occasion later in this report to make use of this obser- 
vation in investigating the effects of a lag in the time required for light 
production from froth particles in the meteor wake. However, first we shall 
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consider the predictions of the model based on the conventional luminosity 
equation. In that case the values of computed in equation ( 13 ) are iden- 

tified as the luminous efficiency factor aefined by equation (l) or ( 3 ) . 

Statistical Procedure 

The iterative procedure for analyzing meteor data described in the 
previous section can be used to compute quantities such as the luminous effi- 
ciency factor, meteoroid density, and mass along the trajectories of individ- 
ual meteors. A machine program has been developed for this purpose and 
applied to 40 of the meteors for which data are given in reference 5 * The 
trail numbers of the 4o meteors are listed in table II. Values of the lumi- 
nous efficiency factor so obtained vary by more than an order of magnitude. 

As discussed in references 5 and 6 , there are a number of sources of random 
error that render the data of individual meteors unreliable. Consequently, 
the authors of those references recommend the use of statistical methods to 
analyze the data of a number of meteors simultaneously. 

In the method adopted here the calculations described in the previous sec- 
tion are carried out for individual meteors and subsequently statistical anal- 
yses are made of the resulting values of quantities such as the luminous 
efficiency factor. Our machine program was expanded to include such processes. 
Tables are formed of the computed and input quantities for all of the meteors 
in the sample. The tables are rearranged in the order of increasing velocity. 
A least squares fit to a straight line is then made of the logarithm of the 
luminous efficiency factor as a function of velocity and the root-mean- square 
deviation computed. This operation is carried out in subintervals of the 
total velocity range since the resulting velocity dependence is not linear 
over the entire range. Figure 3 is an example of the results from such a 
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Figure 3« - Luminous efficiency factor versus velocity. 
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calculation. The circles are the values of the luminous efficiency factor 
from the least squares fit at the mean velocity in the subinterval of velocity 
covered in each case. The vertical bars indicate the root -mean- square devia- 
tion. Each circle in figure 3 is computed from 20 trajectory points from 
about 10 meteors (more or less, depending on the number of trajectory points 
per meteor). It should be mentioned that each circle in figure 3 does not rep- 
resent an independent set of data. Adjacent circles represent two groups of 
data that have 15 trajectory points in common out of their totals of 20 points 
each. Thus only every fourth circle represents independent data. The overlap- 
ping was employed to obtain a better indication of the velocity dependence 
from a limited amount of data. 

In other operations, the tables of computed and input quantities are 
arranged in the order of increasing mass or impact pressure instead of veloc- 
ity. Least squares fits are then made in subintervals of the ordered variable 
to construct plots of luminous efficiency factor versus mass and versus impact 
pressure as shown in figures 4 and 5- Least squares fits and plots of other 
quantities such as mean meteoroid density are also made as shown in subsequent 
figures . 

A 40-meteor sample was selected from among those available in reference 5 
for the foregoing operations. The selection was made on the basis of the 
longest trails and a wide velocity range. Several meteors for which the pro- 
gram indicated densities greater than iron or very large luminous efficiencies 
were replaced by others. The trail numbers of the final 40-meteor sample are 
listed in table II. 

It is stated in reference 5 that the 4l3 meteors for which reduced data 
are given do not represent a random sample of the incoming meteoroid flux. 

For example, the relative numbers in a given magnitude range drop below that 
of larger samples in the faint meteor range. The corresponding distribution 
of our 40-meteor sample is similar to that of the 4l3 meteors of reference 5, 
although it is more erratic because of the small number involved. The bias 
against faint meteors is amplified also because long trails were selected. 

The photographic magnitude range encompassed in our sample is +1.8 to -3.6. 

The computed preentry masses lie between 0.0828 and 3-30 grams. For spheres 
of density 2.8 gm/cm 3 , the corresponding radii are 0.19 to O .65 cm. The com- 
puted radius of the hemispherical froth shell developed during flight in the 
atmosphere averages about 0.7 cm. 


RESULTS AND DISCUSSION 


Results From Statistical Procedure 

As mentioned in the previous section, the results from an analysis of 
super-Schmidt meteor data are presented in figure 3- The scale on the left 
indicates the units used in references 5 and 6 for the luminous efficiency 
factor T 0 p • The scale on the right can be used to compute the equivalent 
dimensionless luminous efficiency T p > which is the fraction of the kinetic 
energy of the vaporized atoms that is converted to light detectable by the 
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blue sensitive film used in the super-Schmidt cameras. This scale was con- 
structed according to reference 6. Additional information on the method of 
conversion is contained in references 1 6 and 17 * Since the quantity read from 
the scale is 10t p /v, it should he multiplied by v/lO to obtain the luminous 
efficiency T p • 

Even though updated material properties have been used, figure 3 is still 
similar to that from the preliminary analysis reported in reference 4 . The 
constant value of T op recommended by Verniani in reference 6 is shown as a 
dashed line in the figure. The values of T op from the present method 
approach Verniani* s value at high velocity, but rise with decreasing velocity 
rather than remaining constant. It was shown in reference 4 that such a veloc- 
ity dependence of Top brings the velocity dependence of the meteoroid flux 
deduced from photographic meteor data into better agreement with overlapping 
radar results. A similar velocity dependence of T op was found for small 
iron spheres in laboratory experiments (ref. l8) . Averages over 10 points at 
a t im e of the 50 data points given in reference l8 are shown in figure 3* The 
necessary conversion of units was based on the assumption that the luminous 
efficiency measured in the experiments was the photographic efficiency T p • 
According to reference 8 the luminous efficiency for stone should be less than 
that for iron by a factor of about 6.5. No readjustment in the results shown 
for iron has been made to convert to an equivalent value for stone. 

Values of T op versus V for three individual meteors are included in 
figure 3 to show typical variations along a trajectory. The individual 
meteors in the present sample cover a relatively small velocity range, but the 
values of T 0 p computed often vary by more than a factor of 2 along a 
trajectory. 

It is known that, due to instrument selection effects and other factors, 
correlations exist between variables such as mass and velocity in meteor data 
chosen at random (ref. 5)- It is therefore of interest to determine whether 
the apparent velocity dependence of Top shown in figure 3 might be basically 
a dependence on mass or other variables. Figure 4 is a plot of T Q p versus 
mass. It appears that the correlation with mass is less than with velocity. 
Typical variations along the trajectory are shown for three individual 
meteors in figure 4 . 

Figure 5 is a plot of T op versus impact pressure, pV 2 , constructed from 
the same meteor data as the previous two figures. It appears that there is 
little or no correlation of T op with this variable for the 40 meteors used 
in the calculation. Typical results for three individual meteors shown in the 
plot indicate that the impact pressure can increase appreciably along a 
trajectory. 

As explained in the Analysis section the present method of meteor analy- 
sis leads to a prediction of the variable mean meteoroid density along a tra- 
jectory. Figure 6 shows the mean meteoroid density as a function of impact 
pressure obtained by averaging over 20 trajectory points at a time. As dis- 
cussed in the section on the statistical procedure , each average is con- 
structed from the data of a number of meteors, depending on the number of 
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Figure 6.- Meteoroid mean density versus impact pressure. 


trajectory points per meteor falling in the range of impact pressure included 
in the average. The computed values of mean meteoroid density versus impact 
pressure are also shown for two individual meteors in figure 6. The results 
for these two meteors indicate that both the mean meteoroid density and impact 
pressure vary appreciably along a trajectory. 

Individual meteors either increase or decrease in density. However, the 
averages over several meteors shown in figure 6 indicate a systematic tendency 
toward higher mean meteoroid densities at the larger impact pressures. This 
indication of a smaller accumulation of froth on the average at higher impact 
pressures can be attributed to the greater likelihood of the froth being blown 
off. It would be difficult to explain this plot on the basis of meteoroids 
not containing solid inner cores. For example, a general fragmentation of 
fragile meteoroids should result in lower apparent densities at higher impact 
pressures, opposite to the tendency shown in figure 6. This observation 
lends credence to the assumption that the solid inner core does not fragment 
for a majority of meteors in the sample. 
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Velocity, V, km /sec 


Figure 7.- Meteoroid mean density as a function of velocity. 

Figure 7 shows the mean meteoroid density as a function of velocity. The 
correlation of density with velocity is less than the correlation with impact 
pressure shown in the previous figure in the sense that the scatter overlaps 
the variation more. The drop-off in mean density at high velocity can he 
related to the lower impact pressures experienced by the higher speed meteors 
on the average. The drop-off in mean density at lower velocity is believed to 
be fortuitous. An unexpectedly large fraction of the low velocity trajectory 
points in the sample occurred at values of the impact pressure less than 
2000 N/m 2 . According to figure 6 , small values of mean meteoroid density are 
to be expected under such conditions. The depressed impact pressure at low 
velocity in our meteor sample is due to a combination of factors. About a 
third of such points correspond to early parts of the trajectories. For sev- 
eral of the meteors involved the later parts of the trajectories were trun- 
cated because they left the field of view of one of the observing cameras. 

The contributions of the larger mean densities that do occur at low velocity 
are suppressed in the logarithmic averages that include the more numerous 
trajectory points showing small mean densities. 
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In reference 6, Verniani reported an analysis of 36 O of the meteors 
contained in reference 5 • It was f ound that the average value of Top v / 
for meteors with aphelia greater than 7 AU was greater by a factor of 2 than 
the average for those with smaller aphelia. The method of reference 6, based 
on the luminosity and dynamic equations, did not permit a determination of the 
mean density Pm alone. It was assumed that the effect was due to a differ- 
ence in density of the two groups of meteoroids. The difference in density 
was justified on the basis of Whipple* s icy comet model. No allowance was 
made for possible changes in mean meteoroid density due to frothing in the 
atmosphere. On this basis T 0 p was found to be independent of velocity if 
the long -period meteoroids were assumed to be less dense by a factor of 1.4 
than the short-period meteoroids. With no difference in density for the two 
groups, it was found that T 0 p would be approximately proportional to V 1/2 . 

In the present method the mean meteoroid density and luminous efficiency 
factor are computed separately. With the aid of information contained in 
reference 19 9 it has been determined that the last three average values of 
mean meteoroid density on the right in figure 7 are largely from meteoroids in 
orbits with aphelia greater than 7 AU. Most of the meteoroids used for the 
averages at lower velocities are from orbits with aphelia less than 7 AU. If 
we ignore the decline in mean density at lower velocity for reasons previously 
discussed, a difference in mean density of about a factor of 2 is found for 
the two groups of meteoroids discussed by Verniani. If this density ratio 
were put into Verniani* s analysis, the luminous efficiency factor determined 
would rise with decreasing velocity rather than remaining constant . As previ- 
ously discussed, the decline in mean density at high velocity can be attrib- 
uted to the greater froth accumulation that is possible at the lower impact 
pressures experienced rather than to a difference in the densities of the 
meteoroids before entry into the atmosphere. In a later section, the effect 
of changes in the assigned parameters is discussed. From that study it can be 
determined that if the density of the solid core were assumed to be smaller 
for long-period meteoroids than those with aphelia less than 7 AU, the 
decrease in mean meteoroid density at high velocity would be accentuated by 
about the same ratio. This would lead to an even faster rise in luminous 
efficiency with decreasing velocity from Verniani *s analysis than would result 
from the density variation shown in figure 7- Likewise, if a lower density of 
the solid meteoroid core were assumed for high-velocity meteors, the luminous 
efficiency factor computed by the method of the present report would drop even 
more than with a constant solid core density (see eq. (l6)). 

A plot of mean meteoroid density versus mass showed little correlation 
and is not reproduced here. The most informative plot of mean density is con- 
sidered to be that showing its dependence on impact pressure in figure 6, 
which was previously discussed. 
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Plots have been made of computed values of the surface temperature at the 
stagnation point versus velocity, mass, and impact pressure. Of these, the 
one indicating the greatest correlation is the plot of temperature versus 
impact pressure shown in figure 8. The stagnation-point surface temperature 
would be a monotonic function of pV 3 except for variations of the aerody- 
namic heat-transfer coefficient, which depends on the velocity, air density, 
and nose radius of the meteoroid. Allowing for evaporation of volatile 
meteoroid components in the aerodynamic energy equation (C2), discussed in 
appendix C, also disturbs this monotonic relationship. The scatter in fig- 
ure 8 reflects these disturbing influences as well as the difference in veloc- 
ity dependence of the impact pressure pV 2 and the airstream energy flux pV 3 . 
Results for two individual meteors plotted in the figure show that the surface 
temperature and impact pressure increase appreciably along a trajectory. 


Sensitivity of Results to Assumptions 

Since the meteoroid properties are uncertain and an idealized meteoroid 
shape has been assumed in the foregoing analysis , it is of interest to deter- 
mine the effect on the computed luminous efficiency factor of changes in 
assigned values of the parameters. For this purpose, it is helpful to derive 
approximate expressions for the dependence of r Q p on the parameters. The 
resulting approximations have been checked by making calculations for a few 
meteors with one parameter at a time varied. A matter of particular interest 
is to determine whether any reasonable change in assumed properties can remove 
the velocity dependence of the luminous efficiency factor indicated by the 
present method. 

By means of approximations of equations (3), (7)* and (8) described in 
appendix D, the proportionality relation 
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has been derived. This expression indicates that if, for example, the 
assigned value of the nose radius parameter A rn is doubled, the computed 
values of r 0 p will be increased by the factor 2 3 / 4 = 1.682. This applies to 
the computed values at each point on the trajectory of an individual meteor as 
well as the averages from the statistical procedure. In figure 9 the above 
prediction is compared with the results of machine calculations for a few 
individual meteors. 

Figures 10 and 11 show similar comparisons with respect to variations of 
the geometrical parameters A rs and Sjp. The shapes corresponding to these 
changes in parameters are discussed in appendix A and are depicted in fig- 
ures 19 and 20. It can be seen from figures 9 to 11 that equation (l6) is a 
reasonably accurate approximation in view of its simplicity. Since the magni- 
tudes of the variations considered a/e in a range that can be expected to 
occur in a sample of meteors, it is clear that the level of the luminous effi- 
ciency factor determined in the previous section is not well determined. How- 
ever, a systematic variation of the parameters with velocity that would remove 
the velocity dependence of the luminous efficiency factor appears unlikely. 

The predictions of equation (l6) with respect to changes in the assigned 
values of the density of the solid core and of the froth density are compared 
with results from machine calculations in figures 12 and 13 . The approxima- 
tion appears to be useful for adjusting the solid core density, but not the 
froth density. For the latter, the variations in the luminous efficiency fac- 
tor are larger than predicted and are in the opposite direction. This failure 
has been traced to the step in the approximation whereby the second term in 
equation (7) is omitted (see appendix D) . This term is ordinarily small com- 
pared to the first term, but increasing by a factor of 3 increases the 

ratio of the two terms by about the same factor so that the effect of the 
second term becomes more important. In addition, the term in question con- 
tains p m as a factor which occasionally becomes large. For this reason, the 
machine calculations also become erratic when a large value is assigned to the 
froth density. It can be argued that large values of this term should be sup- 
pressed when they lead to increased scatter in the computed values of the 
luminous efficiency factor. The subject will not be pursued further here- 
However, it should be noted that if values of the froth density above about 
200 kg/m 3 are of interest, complications not indicated by the approximate 
equation (l6) should be considered. 

Approximations leading to predictions of the effects of changes in other 
assigned parameters are described in appendix D. The proportionality relation 
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has been derived. In figure l4 
results from machine calculations are 
compared with the predictions of equa 
tion (l7) for changes by a factor of 
3 in the assigned value of the ther- 
mal conductivity k- Since the quan- 
tities Cy, k, and C p appear only in 
the combination above in the ablation 
model equations, figure l4 also indi- 
cates the validity of the approxima- 
tion with respect to variations in 
their assigned values. The tangen- 
tial force coefficient C p depends 
on computed quantities such as the 
meteoroid nose radius as well as the 
air density and meteor velocity. 


Gjp is not an assigned 


Although 

parameter, equation (17) indicates 

the effect on T op to be expected when Cy is modified by a constant factor 

such as the constant C in 


due to the reassignment of a parameter 
equation ( C3 ) • 


The equilibrium vapor pressure p v likewise is not an assigned parameter 
However, equation (l7) indicates the effect on T op that would result from a 
change in the level of the computed vapor pressure due to reassignment of a 
parameter in the equilibrium vapor pressure law (eq. (Oil)). Figure 15 shows 
a comparison of the prediction with results from machine calculations for a 
few individual meteors. It is seen that changing the level of the equilibrium 
vapor pressure law by a factor of 10 changes the computed values of T op by a 
factor of about 2. 5- 
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The quantity in equa- 

tion (17) is the heat per unit mass 
required to heat and melt the solid 
meteoroid core; Ty is the tempera- 
ture of fusion. It has been deter- 
mined from machine calculations that 
the effects on T op of changing the 
assigned values of these parameters 
indicated by equation (l7) are 
realistic . 

From information given in appen- 
dix C, it is found that at a typical 
temperature of 2400° K, ordinary 
chondrites have a viscosity greater 
than that of carbonaceous chondrites 
by a factor of 3*4 and a smaller 
equilibrium vapor pressure by a fac- 
tor of 3*6. According to equa- 
tion (17) the net effect of these 
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two differences is an increase in the predicted values of r 0 p by a factor 
of about 1 . 5 . A somewhat smaller effect was found in the machine 
calculations - 

Additional information on the effects of changes in the assigned param- 
eters is given in appendix D. Although the level of the values of r Q p i n 
figure 3 may be in error, no simple reassignment of the parameters has been 
found that would remove the velocity dependence • 


Effect of Froth Ablation in Wake 


The possibility that the mass determined from the luminosity equation 
may differ from the dynamic mass has long been recognized. In reference 
it is observed that if, as a result of aerodynamic forces or thermal stress, 
material is jettisoned from the body in the solid state or if it flows off 

the body in the liquid state, there will be a time lag before the debris 

vaporizes in the wake- The luminosity results from collisions of air and 
vapor molecules and hence does not appear until the material from the body 
vaporizes. If a large fraction of the mass is ablated in the solid or liquid 

form and if the debris follows the parent meteoroid closely enough, the major 

effect may be simply a lag in the luminous intensity relative to the time of 
ablation from the parent meteoroid. The mass determined from the luminosity 
equation would then correspond to the dynamic mass at an earlier time » We 
shall investigate the consequences of such a mechanism on the computed values 
of the luminous efficiency factor. 


The luminosity equation in the usual form with lag effects neglected can 
be written 


m = 



( 18 ) 


In the development described in the analysis section, this equation is used 
to determine T op > the mass and mean meteoroid density being determined from 
the dynamic and ablation equations. Therefore, if equation (l 8 ) is modified, 
the mass and density computed along a trajectory will not be affected. It is 
convenient to retain the definition of T op afforded by equation (l 8 ) while 
defining a corrected luminous efficiency factor ^opc the elation 

/~oo ip r °° ip 

m = 2(1 - B) / *-= dt + 2B / £-? dt (19) 

T opc v 3 J t+e T opc v 3 

The quantity B is a number between 0 and 1 representing the fraction of the 

mass ablated as froth, which is vaporized in a time 9 after ejection from 

the parent meteoroid. The fraction (l - B) of the mass is assumed to be 
vaporized in a time small compared to 9 . It is assumed that a single value 

of 9 and a single value of B in place of an integral over all possible 

delay times can approximate the effect of the delayed vaporization. If so, 
T opc defined by equation (l 9 ) may more closely represent the efficiency of 
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converting the kinetic energy of the vapor molecules to phot ograp hah le light 
and may be less variable than T 0 p • 


From equations (l8) and (l9) the mass before entry into the atmosphere is 
found to be 


“l_oo = 2 




T V" 

-°° 1 op 


dt = 2 




T V*“ 
' opc v 


dt 


If T 0 p C is less variable than T 0 p, its use in the relation 



to determine the meteoroid flux in space is preferable to use of T 0 p in the 
same equation. Determinations of T 0 p C from meteor data may also snow less 
scatter than T 0 p it part of the variability of T 0 p is due to an intensity 
lag. It -will be shown that for the 40 meteors analyzed in this report the 
scatter in T 0 pc l ess than that in T 0 p • 

In appendix E an approximate relationship between T 0 p C and T 0 p is 
derived in the form 

T opc = T op (^- + 

Since m/m = d2n(m)/dt and T 0 p are computed along the trajectories of indi- 
vidual meteors in the machine program described in the analysis section, this 
expression can be used to compute values of Tope ^ ike product B Q can be 
evaluated. We shall first estimate B 0 from a consideration of the behavior 
of froth particles in the wake and later describe a method for deducing the 
value from the meteor data. 

From the calculations for individual meteors previously discussed, it is 
found that most of the material ablated from the front face of the solid inner 
core is removed by flowing to cooler regions rather than by vaporization. 

This behavior is discussed further in appendix C. The mode of ablation of the 
froth produced from this flow is not determined by the ablation model 
described in the analysis section. However, since the froth would be in a 
cooler region than the front face of the solid core, it can be assumed that 
most of the froth is ablated in the solid or liquid form rather than by evapo- 
ration while attached to the parent meteoroid. For that reason we assume that 
the quantity B in equations (19) and (20) is approximately equal to 1. 

The solidified chunks of froth in the wake would also be expected to 
ablate by liquid runoff with the drops subsequently breaking up as a result of 
aerodynamic pressure. When the melted particles produced in this process are 
small enough to avoid disruption by virtue of surface tension, the subsequent 
reduction in mass can occur only by vaporization. When the heat sink repre- 
sented by a progressively exposed cooler interior is absent, the surface 
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temperature is elevated and rapid vaporization occurs. Our calculations 
indicate that most of the meteoroid material reaches the vaporized state by 
this chain of events rather than by direct vaporization from the parent 
meteoroid (see appendix C) . 

An approximate expression for the time required for a solid or liquid 
particle to ablate most of its mass is derived in appendix E in the form 

( 21 ) 

where 


9 = 


6£p; 


2/3 1/ 3 

mf fo 




£ specific heat of ablation 

Pjnf particle density 

m fo initial mass 

Sjn shape factor 

Cjj heat-transfer c oef f ic ient 

pV 2 impact pressure 

V particle velocity 

All quantities in this expression except the initial mass m^ Q can be readily 
estimated. We are concerned with a condition that Opik (ref. 3) called an 
abnormal environment in which the particles involved are injected by ablation 
from a larger body at an altitude lower than they would attain if injected 
directly from outside the atmosphere. The impact pressure and velocity do not 
change appreciably during ablation under this condition. The values 
pV 2 = 3400 N/m 2 and V = 2X10 4 m/sec corresponding to an early part of the 
trajectory of meteor 4464 of reference 5 will be used to estimate 0. 


We shall first consider the liquid drops held together by surface ten- 
sion. From a formula in reference 3 it is shown in appendix E that for the 
above conditions the radius of such particles , considered to be spheres, is 

Ti = 5-1X10 4 meter 

With the froth density p^-p = 100 kg/m 3 , the corresponding mass is 

m z = 5-5X10" 8 kg 

Substituting this value for m^ Q and £ = 8.51x10 s , p mf » = 100, Cjj = 1, 
pV 3 = 6.8X10 7 , S f = 1.208 in equation (2l) yields 

9^ = 0.051 sec 

as the time required to vaporize such liquid drops. 
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The initial mass of solidified froth particles can he estimated if it is 
assumed that the sloughing of such particles is observable as fluctuations in 
the light curve. From light curves in reference 20 and the ablation rate of 
meteor 4464, a value of 


m^ = 3X10 " 5 kg 

has been arrived at (appendix E) . The radius of a spherical particle of this 
mass and density = 100 kg/m 3 is 

r f = 4. 2X10 " 3 meter 

Since liquid ablation rather than vaporization would occur , the specific heat 
of ablation is taken to be l = 1.884x10 s ( from table I) . With other 
constants evaluated as before, the resulting lag time from equation ( 2 l) is 

9^ = 0.092 sec 

Since the light is not produced until the ablated liquid drops vaporize, the 
total lag time is 

0 = 0 _p + ©2 = 0*1^3 sec 


The froth particles are smaller and less dense than the parent meteoroid. 
They would therefore lag behind and produce most of the light some distance 
back in the wake . In appendix E an expression for the lag distance is 
derived in the form 


AL = 


1 H 

2 C H V 


0 


( 22 ) 


The displacements of the liquid drops and froth particles according to this 
relation with Cp = 2, C H = 1, V = 2X10 4 meters are 

AL^ = O.O 65 km (liquid drop) 

ALp = 0.026 km (froth particle) 


The total displacement of the region of light production from the parent 
meteoroid is approximately the sum of these 

AL = 0.091 km (total) 


For comparison, the distance covered by a typical meteor between camera shut- 
ter breaks is 0.3 km, which is 3 to 10 times larger than the above predicted 
lag distances. 

It is interesting to determine whether there is evidence in the meteor 
data of an intensity lag time of the order of 0.1 second as predicted in the 
foregoing discussion. Toward that end, equation (20) can be rearranged in the 
form 
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1 


1 


B0 m 


(23) 


op 


opc 


T opc 111 


The quantities T 0 p and m/m are computed at each trajectory point in the 
method described in the analysis section. If T Q p C and B0 /t opc are assumed 


to be constant along the trajectory, the computed values of T op and m/m at 

two trajectory points can be used in 
equation ( 23 ) to solve simultaneously 
for t op c and B0/Tqp C . By this means 
values of B0 as well as T 0 p C can 
be computed from the meteor data. 
Figure 16 shows values of B 0 thus 
determined for points on the trajec- 
tories of 15 meteors. The computed 
delay times are plotted versus 
10 8 /pV 3 to check the dependence on 
this variable indicated by equa- 
tion (2l) . A majority of the points 
fall reasonably close to the esti- 
mated value 0 = 0 . 1^3 sec arrived at 

Figure 16 .- Comparison of intensity lag time ln the preceding discussion. (The 

from meteor data with estimates (B = 1 ). estimated value Of B is 1.0.) 



In appendix E approximate formulas for the delay time as a function of 
10 8 /pv 3 are derived to obtain 

e f -= o.o 63 xio 8 /pv 3 (24) 

This is the time required to liquefy a froth particle of initial mass 
nifo = 3X10 " 5 kg (shown as a dashed line in fig. 16) . The time required to 
vaporize liquid drops is approximately (appendix E) 

6 Z = 0.0235(lO®/pV 3 ) 2 (25) 

shown as a broken line in figure 16. The total time required to liquefy the 
froth particles and vaporize the resulting drops is approximately Of + Oi 

shown as a full line curve in figure 16 . Since the initial mass m ^ 0 may 

vary considerably from the assumed value of 3 X 10“ 5 kg, the comparison of the 
computed points with the full line curve is satisfactory except for several 
low points. The low values of B 0 could be explained by an unusually high 
volatility of the meteroid material such that B is small compared to 1 for 
these meteors. In that case, vaporization may dominate over liquid runoff for 
the froth particles. The tendency of values of B 0 computed from the meteor 
data to lie above the full line curve at small values of 10 s / pV 3 may be due to 
the presence of larger froth particles' or a greater froth density than 

assumed, or a combination of both. Since figure 16 shows that B0 varies 

from meteor to meteor, it is considered preferable to determine its value for 
each meteor from the data for purposes of evaluating T 0 pc (rather than 
utilizing a single value of B 0 for all meteors) . 
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As previously mentioned , our machine program computes values of T 0 p and 
m only at interior points on the trajectory of each meteor. Since calcula- 
tions at two or more trajectory points are required to evaluate T 0 pc anc ^ 
from equation ( 23 ) > such evaluation can be carried out only for meteors with 
data given at 4 or more trajectory points in reference 5* For this reason, 
T 0 pc can computed for only 15 of the 40 meteors listed in table II (those 
marked with the letter a) . The statistical procedure previously discussed 
has been used to determine the dependence of T 0 p C on velocity. When Tope 
could not be computed, it was set equal to T 0 p • The results are shown in 
figure 17. Comparison with figure 3 indicates that the scatter in computed 
values of T 0 p C is less than the scatter for T 0 p, although the velocity 
dependence of the two is similar. There is a small amount of flattening of 
the velocity dependence at high velocity and a general lowering of values at 
all velocities. 
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Figure 17 . - Luminous efficiency factor with intensity lag correction. 
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The results for iron from laboratory tests (ref. l8) are also shown in 
figure 17* The prediction of reference 8 that the luminous efficiency of 
stone should be less than that of iron by a factor of about 6.5 is more nearly 
borne out by these results than the uncorrected values shown in figure 3* The 
remaining discrepancy may be due to a failure of the assumption that the 
photoelectrically measured luminous efficiency of the tests is equal to the 
photographic luminous efficiency. 

In reference 6 it is noted that calculations of T op based on an 
assumed constant density show an apparent increase of Top along the 
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trajectory for many individual meteors. The effect is ascribed to progressive 
fragmentation of the virgin meteoroid material. Verniani concludes that, con- 
sequently, only the values of T 0 p from the early parts of the trajectories 
are reliable. Increases in T 0 p along the trajectory also result from the 
present method for many individual meteors. In contrast, the values of T 0 pc 
are nearly constant along the trajectories in all cases. Consequently, appar- 
ent increases of T 0 p in references 5 a nd 6 may be due, in part, to the inten- 
sity lag effect. Both decreases and increases of mean meteoroid density along 
the trajectories are noted from the present method. When decreases in mean 
density occur, the evaluation of T 0 p "by the method of reference 6 (substitu- 
tion of a constant value of % in the ratio T 0 p/p^) would also lead to 
apparent increases of Top • 

CONCLUDING REMARKS 


An ablation model for solid stone meteoroids allowing for a decrease in 
mean density due to frothing has been used to analyze the data of 4o meteors. 
Material properties appropriate for carbonaceous chondrites were assigned. 

The photographic magnitudes of the sample ranged from +1.8 to -3*6. The pre- 
entry masses computed were between 8xi0~ 5 and 3-3X10” 3 kg, corresponding to 
radii between 1.9X10" 3 and 6.5X10 -3 meter for solid stone spheres. The com- 
puted mean densities during atmospheric flight ranged from 300 to 1000 kg/m 3 , 
corresponding to flanges of froth as thick as 7X10- 3 meter. 

The luminous efficiency factor was computed along the trajectory of each 
meteor and averaged over several meteors in small velocity intervals. The 
resulting averages agreed with the constant value recommended by Verniani at 
velocities above 35 km/sec. For velocities below 25 km/sec, T 0 p was found to 
rise to about six times Verniani f s value. 

The effect of a delay in light production resulting from the time 
required to vaporize chunks of froth in the wake has been investigated. The 
scatter and the level of computed values of the luminous efficiency factor 
were reduced slightly by correcting for this delay. The effect can account, 
in part, for apparent increases of T 0 p along the trajectories previously 
attributed to fragmentation. Such increases can also result from the assurp- 
tion of a constant mean density in cases where the computed mean density 
decreases due to accumulation of froth. It was found that the data of the 40 
meteors analyzed can be explained on the basis of originally solid meteoroids 
that do not fragment except by sloughing froth produced during the flight. 

The properties such as viscosity and equilibrium vapor pressure assigned 
were chosen to correspond to carbonaceous chondrites. The density of the 
solid inner core was assumed to be 2800 kg/m 3 and the density of the froth 
produced during flight was taken to be 100 kg/m 3 . An investigation of the 
effect of changes in the assigned material properties was made. Uncertainties 
in the properties could alter the level of the computed luminous efficiency 
factor but not its velocity dependence. Assignment of properties appropriate 
for ordinary chondrites would not affect the results significantly except for 
an increase in the predicted stagnation-point surface temperatures. However, 
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many of the above conclusions would appear to break down if a froth density 
more than twice the assumed value of 100 kg/m 3 were assigned. 


Ames Research Center 

Rational Aeronautics and Space Administration 

Moffett Field, Calif., 9^035 > June 5 , 1968 
129 - 01 - 02 - 09 - 00-21 
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APPENDIX A 


DISCUSSION OF ASSUMPTIONS 


In addition to the assumptions listed in the analysis section, specific 
values must be assigned to geometric parameters and constants in the material 
property laws to carry out an analysis of meteor data. For example, the first 
assumption can be stated in the form 

r n/ r = ^rn (Al) 

For most of the calculations, A rn = 1 has been used corresponding to a hemi- 
spherical front face. However, the effect on the results of changing the 
ratio to a larger value corresponding to a blunter shape has been investigated. 

The second assumption can be expressed as 

r s /r = A rs (Q s /Q) l/3 (A2) 

where A rs is a constant, Q s is the volume of the solid core, and Q, the 
total meteoroid volume. The origin and meaning of this relation requires fur- 
ther explanation. Such a rule would apply, for example, if the solid and 
total volumes were both spheres or hemispheres. Equation (A2) was originally 
derived for the case of a hemispherical overall shape and a cylindrical solid 
volume with the length of the cylinder equal to its radius. In the limit of 
very small r s /r it can be shown that A rs = ( 2 / 3) 1/3 for that case. This is 
the value used in the calculations, but the effect on the results of changes 
in the value of A rs has been investigated. When there is no froth, both 
r s /r and Q s /Q should equal 1, whereas the above value of A rs leads to 
r s /r = 0.874 at Q s /Q = 1. This 13 -percent discrepancy can be used as a basis 
for judging the error that might arise from the use of assumption 2. 

The third assumption is imposed largely for reasons of mathematical 
expediency and leads to simplification in the analysis of the viscous flow of 
the thin layer of melted material covering the solid. The presence of meteor- 
oid rotation or oscillation about an axis not alined with the flight direction 
is denied by this assumption. Such motion might be expected for about one- 
third of a random sample of meteoroids; however, the amount of scatter in the 
results due to this assumption is unknown. 

The fourth assumption is a statement of the modes of ablation of the 
solid core that are considered in the present investigation. The last assump- 
tion, regarding the composition of the meteoroid, supplies necessary input for 
the calculations. Most of the results are based on material properties appro- 
priate for carbonaceous chondrites. However, the effect on the calculations 
of variations in properties has been investigated. It has been found that the 
values of luminous efficiency factor and meteoroid mean density computed are 
essentially unchanged when properties appropriate for ordinary chondrites are 
assumed, although higher surface temperatures result. 
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The implications of the geometrical assumptions should perhaps he further 
illustrated. For most of the meteor data to which our ablation model has been 
applied, the value of the computed ratio of radii r s /r lies between 0.3 and 
0 . 8 . It is of interest to determine the changes in shape that correspond to 
variations of this ratio as well as changes in the shape factor Sp defined 
in equation (4) . For that purpose it is expedient to make additional geomet- 
ric assunptions . However, it is worth noting that these additional assump- 
tions are made only for purposes of illustration and are not necessary for 
application of the ablation model to meteor data . 


Since the frontal area is assumed to be circular, equation (4) can be 
rearranged in the form 



or 



(A3) 


where Sp^ is the shape factor for a hemisphere 




Figure 18. - Sketch of circular cylinders with 
spherical ends used in alternative interpre- 
tation of model geometrical relations. 


The construction of shapes corre- 
sponding to various shape factors and 
ratios of radii r s /r is greatly sim- 
plified with the assumptions of axial 
symmetry and spherical ends on the 
solid and froth cylinders as illus- 
trated in figure l 8 . In that case the 
volume of the solid can he written 

Qs = rtr i J s 

Combining this with equation (A 3 ) and 
rearranging yields 

2s = 2 / s FhV/ 2 Qs/Q 
r 3 V S F / (r s /r ) 2 

Similarly , the total volume is 
given by 

Q = Jtr 2 Z f - Jtr|(Z f - Z s ) 
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Combination with equation (A 3 ) and substitution of equation (A4) leads to 


ll 

r 


2 

3 


%h \ 3/2 

ftr ) 





(A5) 


The relation arising from assumption 2 given in equation (A 2 ) can be 
rearranged in the form 

Qs.l 

Q 2 


•r s \s 


(A 6 ) 


S F = 3.05 S F = 1.919 


S F = 1.208 



(The value A rs = ( 2 / 3) 3 used in 
the analysis of meteor data has been 
inserted.) Substitution of this into 
equations (A4) and (A 5 ) leads to the 
relations 



3/ 2 


£s 

r 


(A7) 


and 


1 - 2 (la 
2 /%hV /2 2 


r 3 \ S p / 


1 ■ 1 r 


'r s ' 2 


(A8) 


Figure 19 shows the axially sym- 
metric shapes with spherical ends 
corresponding to various shape fac- 
tors Sp and various ratios of solid 
volume to total volume Q s /Q. The 
sketches are generated from equa- 
tions (A6) to (A8) with the added con- 
dition, nose radius r n equal to 
base radius r . The column of three 
sketches on the right corresponds to 
a shape factor Sp = 1.208, which is 
the same as that of a sphere . The 
shape factor Sp = 1.919 of the mid- 
dle column is the same as that of a 

hemisphere, while the value Sp = 3-05 of the left column corresponds to what 
might be termed a thin shell. The three sketches in the top row correspond 
to a ratio of solid to total volume Q s /Q = 0.0405. The middle and bottom 
rows correspond to Q s /Q = 0.1875 and 0-768, respectively. Most of the 
meteors analyzed were found to fall in the range of Q s /Q covered here. The 


Figure 19.- Sketches of axially symmetric shapes 
corresponding to various values of shape fac- 
tor Sp and ratio of solid to total volume 
Q s /Q when the front face is a hemisphere. 
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dotted lines in the middle row of sketches show possible rearrangements of the 
rear shapes that are allowed within the set of assumptions for the ablation 
model previously listed. As far as the geometrical assumptions of the model 
are concerned, even nonaxi symmetric rearrangements of solid and froth volumes 
are allowed as long as the relative values of the volumes are not changed and 
the front faces are not altered. 


Q s 

= - 0405 



r n 


Other allowed shape changes are illustrated by the fact that for the mid- 
dle column the froth and solid volumes could be packed within a hemisphere 
without altering the results from analysis of meteor data. Similarly, the 
volumes in the right column could be fit into spheres without altering results 

based on a spherical shape factor. 
The middle column corresponds to 
the value of shape factor used in 
most of the calculations in this 
report, although the effect of a 
change in shape factor is investi- 
gated. The ratio of volumes 
q s /q is determined at each point 
on the trajectory of each meteor 
from the meteor data with the aid 
of the assumptions of the model. 


fp 

s s 


r s 


= .5 


Qs 

Q 


= .1875 


a 1 

s 

\Zi 





Figure 20.- Axially symmetric shapes corresponding 
to various degrees of nose bluntness and ratios 
of solid to total volume (Sp = 1 * 919 ) * 


Figure 20 illustrates the 
effect on the geometry of changes 
in the ratio of nose radius to 
the radius of the frontal area 
r n /r. Again, rearrangements of 
the volumes are allowed that pre- 
serve the front faces as shown. 
The right column with the hemi- 
spherical front face corresponds 
to the ratio r n /r = 1 used in 
most of the calculations, but the 
effect of varying the ratio is 
investigated . 

Next we shall derive or list 
the major mathematical relations 
generated from the ablation model. 
Auxiliary equations arising from 
the aerodynamic and material 
property laws are given in 
appendixes B and C. 


Substituting Q = m/pj^ 


into equation 


(A3) leads to 


r = 

Equation (Al) can be written 


^_y / 3 / ^y/ a 

s2V \%V 


r n = A rn r 



(A9) 

(A10) 
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If values of the shape factor Sp and geometric parameter A rn are assigned 
and the meteoroid mass m and mean density are known, the radius of the 
frontal area r and nose radius r n can be computed from these equations. 


The relation 

PmQ- - PmsQs + Pmf (Q " Q s ) 

expresses the fact that the total meteoroid mass is equal to the sum of the 
masses of the solid and froth parts. Rearrangement leads to 

Qs = p m ~ p mf 
^ p ms ” p mf 

and substituting this into equation (A2) yields 


r s 

— = A. 


rs 


m 


p mfV /3 


ms 


p mf. 


(All) 


If values of the froth density p^, solid density P ms , and geometric param- 
eter A rs are assigned, the radius r s of the frontal area of the solid can 
be computed when r and p^ are known. 


The formula 


m 

p m 


m s m - m s 
p ms p mf 


arises from the equality of the total meteoroid volume to the sum of the 
volumes of the solid and froth parts. Differentiation and rearrangement leads 
to the expression 


ms f p m " Pmf\ . m PmsPmf 

( m + 

Pm Vms - Pmf/ (p ms - Pmf )p 3 


(A12) 


which relates the rate of melting of the solid mass to values and derivatives 
of m and p^ 

Combining equations (A$?) and (All) with the value = (9^/4) 3 

inserted leads to the equation 


A s = 



PmfV /3 

P mf/ 


(A13) 


for the frontal area of the solid core . The rate of melting of the solid core 
per unit area and time is given by 


co 


s - 



(Al4) 


Substituting equations (A12) and (A 13 ) into this and rearranging leads to 
equation (7) of the text. 
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For use in appendix B the relation 



is derived by combination of equations (A 9 ) and (A10) with the value 
Sjph = (9*A) 1/3 inserted. 


4o 


(A15) 



APPENDIX B 


VISCOUS FLOW EQUATIONS 


The relations governing the viscous flow of the melted material covering 
the front face of the solid are presented here in a form applicable to analy- 
sis of arc- jet tests as well as meteors. The configuration under 
consideration is shown in figure 2. 


Basic Relations 

An expression for the conservation of mass with respect to an annular 
control surface at the shoulder of the solid can be written 

-p ms y = 2Pm ^- c Si(T) + S(T) (Bl) 

where 

y rate of recession of the solid surface at the stagnation point 

T surface temperature at the stagnation point 

u c tangential surface velocity of the melt at the control surface 

6 C thickness of the melted layer at the control surface 

gl (T) factor to allow for (l) nonlinearity of the velocity profile at the 
control surface, (2) nonuniformity of the recession rate, and 
(3) nonuniformity of the surface temperature 

£>(t) rate of evaporation per unit area 

The structure of equation (Bl) can be understood if g-^T) is set equal to 1. 
Then the left side represents the rate of production of melted material per 
unit area, considered uniform on the entire front face covering the solid. 

The first term on the right represents the rate of flow of froth through the 
annular control surface at the shoulder of the solid divided by the frontal 
area of the solid, Jtrf . The last term represents the rate of evaporation per 
unit area, considered uniform over the frontal area. 

A question arises as to whether the mass flow through the control surface 
(the first term on the right) should be based on the density of the solid p ms 
or the froth density Pmf In the arc- jet tests, the condition of the models 
after the tests indicated the material to be frothed in this region. However, 
this may be due to the sudden drop in pressure resulting from turning off the 
airstream at the end of the test. In checking equation (Bl) against the arc- 
jet tests, this uncertainty is not important because the combination Pmf^c 
is a measure of mass per unit area which can be determined unambiguously 
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whether the frothing occurred during or after the test. Also, in application 
of the equation to meteor data, it will be seen that Pmf and 5 C are elimi- 
nated by combination with the other conservation equations. Therefore, the 
value of the density assigned for the flow in this region is immaterial as 
long as the same value is used consistently in the other conservation 
relations. 


The momentum equation at the control surface at the shoulder of the solid 
can be written 

H f (T)u c g a 


= | C F pV 


■2 

r n 


(B2) 


where 


P^(t) viscosity of the froth evaluated at the stagnation-point temperature 

g 2 factor to allow for the nonlinearity of the velocity profile and 
nonuniformity of the temperature 


Cp skin-friction coefficient 


The viscosity of the froth is used here to be consistent with the use of the 
froth density in the previous equation. A value of g 2 = 1 would correspond 
to a linear velocity profile. A value of Cp - 2 would correspond to free 
molecule flow, but departures from this condition occur in the analysis of the 
meteor data in reference 5* The application of equation (B2) to the results 
of the arc-jet tests reported in reference l4 led to a prediction of values of 
viscosity consistent with independent measurements. Additional unpublished 
measurements beyond those reported in reference l4 were made by Charles 
Shepard at Ames Research Center for this purpose. 


An energy equation can be written for a steady state with the assumption 
of uniform conditions over the frontal area of the solid in the form 


k(T - T f ) 

-p ms y£fg 3 = — 5^ 

where 


(B3) 


heat per unit mass required to heat and melt the virgin material 
g 3 factor to allow for correction of nonuniformities and nonlinearities 
k thermal conductivity 
Tp temperature of fusion 

Up to the present, the correction factor g 3 has been set equal to 1 for 
analysis of meteor data. Equation (B3) predicts surface temperatures in the 
arc-jet tests ranging up to about 2000° K. The surface temperatures were not 
measured with sufficient precision to check these predictions accurately. It 
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has been checked that a steady-state equation of the above type is a 
reasonable approximation for the visible part of meteor trajectories if the 
meteoroids involved are stone . 


For application to meteor data, equations (Bl) to (B 3 ) are combined to 
eliminate u c and & c in the form 


where 


\ V p mf 

ft) Mf(T) 


(T - T f ) 2 G = - 


( Pms^ ) ( Pms^ r n 

CppV 2 


G = 


§i 

S2&3 


(b 4 ) 

(B5) 


Available information on the properties of stone at elevated temperatures 
applies to the unfrothed liquid state. We have assumed that the thermal con- 
ductivity of the froth is approximately the same as that of the liquid because 
of the conductivity of vapor in the bubbles. The viscosity of the froth is 
expected to be lower than that of the liquid, and we have assumed 


p f (T) ^ p(T) 
p mf p ms 


(B 6) 


Equation (8) of the analysis section follows from equation (B 4 ) by substitu- 
tion of equations (AI5) and (b 6) and the relation 


pm S y = - w s 


(B 7 ) 


Correction Factor for Nonlinear Velocity Profile 

The contribution of the nonlinearity of the velocity profile to the 
correction factor can be determined by setting 


and 


u c 5 cSi = / ^ dy 

J o 


^ C |X f( T )g2 du 

6 C Pf (T) dy 


(B8) 


(B9) 


where u(y) and T(y) are the fluid velocity and temperature at any point y 
in the profile. If replacements according to these relations were made in 
equations (Bl) and (B 2 ) , they would then be exact with respect to nonlinearity 
of the profile. By partial integration and rearrangement, equation (b 8) can 
be written 


u c 5 c§ 




du 

dy 


(& c ~ y)dy 


^3 



If equation (B9) is substituted into this, the relation 


Si i r 5 c 


§2 85 


Mf(ff) 


(8 - y)dy 


can be derived. 


In accord with reference 21 the viscosity law is taken to be of the form 


Mf(T) n(T) i/T 

— — o M e 


Pmf Pn 


The temperature profile is assumed to be linear so that 


T = T - (T - T f )^l - jr-) 


(Bll) 


(B12) 


Substituting the last two relations into equations (BIO) and (B 5 ) with g 3 = 1 
yields 

G(T) = f X exp-f- — j T) dr) (B13) 


G(T) = / expj T| dr) 

J o 4 T[T - (T - T f )rj]-J 


r] = 1 - (y/s c ) 


In the machine program, G is evaluated by 10-point Gaussian quadrature with 
the interval of integration divided into two parts. 


Correction to Account for the Effects of a Pressure Gradient 
and Meteoroid Acceleration 

In the present notation, equation ( 9 ) of reference 12 can be written 

(Bl4) 


(B15) 

The first term in P 11 is associated with the augmentation of the viscous 
flow due to a pressure gradient and the second term provides an estimate of 
the reduction in mass removal rate due to the meteoroid acceleration. 

Strictly, equation (Bl4) applies only at the stagnation point, but to estimate 
the effects under investigation we shall apply it up to the control surface at 
the shoulder of the solid inner core. Integration of equation (Bl4) from 
x = 0 to the control surface at x & r s yields 


5u _ <*t w f 6 + pH f 5 y d y 

S dx Jo h Jo P 

where 

p ms r n C lA \ 2pV 2 

4m J r 2 
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I 


U C - T wc 




ydy 


(B16) 


Equation (B2) follows from this equation if the second term is neglected, p. 
is considered constant in the integral over y , and with the substitution 


1 wc 


I c F pv 2 


r 


s 


r 


n 


The same operation applied to both terms of equation (Bl6) , with m/A 
replaced by (2/3)Pm r in equation (B15 ) ? leads to 


^ f (T)u c g 2 ± 


= ± c-ppv 2 ia 


1+1 


- n 


1 p p ms r n 

“ 8 C D 

H m 


25 c ' 
C F r n- 


(B17) 


(Bl8) 


Comparison of equations (B2) and (Bl8) shows that the correction under 
consideration can be imposed if Cp is replaced with 


Cp _ = Ctji 

^eff i 


1+11 

\ 


a 


1 p p ms r n 

8 D V 


28 r 


C F r n J 


(B19) 


Values of £> c computed from the meteor data with the aid of equations (B 3 ) 
and (B7) are less than 10" 4 meter. Also^ Cp > 0.2, Pms/Pm < 9> 
r n > 4x10“ 3 meter for the meteors analyzed. With these extreme values the 
correction from the pressure term alone is 2& c /(Cpr n ) < l/4. When the accel- 
eration term is included, the correction term is - 19 / 32 , which is numerically 
larger. In the latter case the quantity in brackets is 13/32 rather than 1.0 
as it is in the absence of pressure gradient and acceleration effects. In 
other words, 

- 12 


CF eff = % °F 


(B20) 


The correction represented by equation (B19) has not been applied in the 
analyses of meteor data of this report. However, an estimate of the error 
from omission of the correction can be obtained from equation (17) simplified 
to 

- 3/4 . 

T op « Cp (B21) 

Thus a change in Cp by the factor 13/32 indicated by equation (B20) would 
increase the computed value of T 0 p by the factor (l3/32 )~ 3/4 = 1-965- Most 
of the computed values of T 0 p would not receive so large a correction since 

this value is extreme. However, some of the scatter in the results may be due 

to omission of the correction. 


Comparison With Another Calculation Method 


The overall accuracy of the foregoing approximate procedure has been 
checked by comparison with results from the more exact methods of reference 12. 
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This was done by constructing artificial meteor data (values of h, V, V, and 
m along a trajectory) using the machine program described in reference 12. 

The machine program of the present report was then used to analyze that data 
with set equal to p ms and other material properties equated. Only one 

step of the iteration described in the analysis section was performed because 
the result of the comparison rests on whether the value of R from equa- 
tion (12) is equal to 1 in the first pass. For the cases considered, the 
values of R computed differed from 1 by less than 20 percent. This indi- 
cates that for a meteor the luminous efficiency computed is expected to be in 
error by less than 10 percent due to the use of the mathematical approxima- 
tions described in this and the next appendix. It is to be noted that the 
approximations lead to a considerable reduction in the machine computing time 
required to analyze the data of large numbers of meteors. In addition, it is 
not known whether a more elaborate ablation model leading to the evaluation of 
the luminous efficiency from meteor data is feasible. 
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APPENDIX C 


AUXILIARY RELATIONS FROM AERODYNAMIC AND PROPERTY LAWS 


Aerodynamic Relations 

The drag coefficient Cj) in equation (6) is evaluated by the formula 

Cp = 1 + e -i* s 37Xio6 _ r (ci) 

taken from reference 12. The air density p and its ratio to the sea level 
value p are evaluated according to reference 11. 

The temperature T at the stagnation point is found by iterative solu- 
tion of the energy balance relation 

\ C H PV 3 = eoT 4 + £ V W(T) + C ms £ f u> s (C2) 


where 

Cjj heat-transfer coefficient 

e surf ace emi s sivity 

£ v heat per unit mass required to heat, melt, and vaporize the solid 

C ms fraction of solid that vaporizes readily or is explosively spalled as 

it reaches the surface 

w(T) rate of vaporization (eq. (CIO)) 

heat per unit mass required to heat and melt the solid 

co s rate of melting of solid core mass per unit area and time (eq. (Al4)) 

The term on the left is the heat input from the air stream. The first term on 
the right represents the heat loss due to surface radiation and the second 
term, the heat utilized in vaporization of the refractory part of the material ■ 
The last term allows for the presence of a volatile component such as water 
that will evaporate readily when it reaches the surface and which has a negli- 
gible latent heat of vaporization after being heated to the melting point of 
the solid. A fraction of the remaining refractory material might also be 
carried off by the boiling action from the volatile components. The calcula- 
tions of this report are based on the properties of carbonaceous chondrites 
which contain as much as 20-percent water. However, a majority of such 
meteorites contain less than 20-percent water and we have adopted a value of 
c ms " 0-1 5- 
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The vaporization term containing a j(t) is dominant in the calculations 
for the meteors of reference 5* Because of the steep temperature dependence 
of that term, the temperature computed from equation (C2) is insensitive to 
the values of all parameters appearing in it except those involved in the 
functional form of co(t) itself. For example, it can be shown that with the 
parameters assigned in table I an increase in C or a decrease in by a 

factor of 2 increases the temperature computed by less than 60° K. Changes in 
the parameter e or C^g in equation (C 2 ) by a factor of 2 would have less 
effect. For the calculations in this report we have set e = 0.6, 

= 8 . 5 IXIO 6 (m/sec ) 2 and = 1.884xi0 6 (m/sec) 2 . The values of all 
parameters are listed in table I. 

The additional terms (l - C ms KfW s - £f w (T) would be required in equa- 
tion (C 2 ) for computing an average temperature over the entire solid frontal 
area. This would account for the fraction of the input energy utilized to 
heat and melt the solid material removed by the viscous flow. The resulting 
decrease in temperature below that computed from equation (C 2 ) is not evalu- 
ated in the present model. However, the effect can be accounted for approxi- 
mately by setting C ms in equation (C2) equal to an artificially large value. 
The insensitivity of the results to an increase of C ms by a factor of 2 
above the assigned value is considered to be an indication that computation of 
the temperature drop is not essential. 


The heat -transfer coefficient Cy is evaluated by formulas taken from 
references 7 and 12 as follows: 


C -51 
°Hcru “ m 


V 


•J P r n 


1 - CN 


1 + 




aV e 

V.eoT 4 + (; v co. 


exp (-7 .2X10 8 /pr) | 


(C 3 ) 


<=Beu * feT' 0 ^- 7 - 2 ” 0 " 8 ^ 


oV n 


c„... = 0.6x10" 3 feY e -i-=xi°- 8 /pi- 


"Huu 


^Hu “ C'Hcvu + ^Heu + ^Huu 


-2 C t 


C H = c Hu e Hu + (1 - e vnu )‘ 


~ c Hu \ 2 


(C4) 

(C5) 
(c 6) 
(C7) 


The units in this report are MKS except where noted. Values of the required 
parameters used in the calculations are listed in table I. The quantity 
C^cvu as a c °nvective heat-transfer coefficient that allows for the presence 
of ablation vapor. It is not necessary to include the radiative heat-transfer 
coefficients Cft eu (equilibrium) and C^u (nonequilibrium) for the analysis of 
meteors corresponding to meteoroids below 1 gram in mass. These terms 
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contribute significantly only in the case of bright meteors. Equation (C7) is 
a bridging formula applicable in the entire range of continuum., transition, 
and free molecule flow. This relation, taken from reference 7 > compares favor- 
ably with a bridging formula derived in reference 12. In the analysis of the 
super-Schmidt meteor data of reference values of Cjj ranging from 0.1 to 
1.0 were computed. 


The rate of vaporization per unit mass for free molecule flow WpM is 
given by the Langmuir equation 


" 0-1383 


s/M p v (T) 


(C8) 


where 


M molecular weight 

p v (T) equilibrium vapor pressure (eq. (Cll)) 


When governed by diffusion in the presence of a viscous boundary layer, the 
rate of vaporization can be approximated by 


co 


d “ 


34 .4mc Hcvu p v (t) 


Pv(T) ~ 

peFT 


which is a modified form of a relation given in reference 12. 


(C9) 


A bridging formula for the general case derived 
relation 


1 = 1 
O) 0)^ 


+ 




in reference 12 is the 

(CIO) 


The equilibrium vapor pressure in equations (c8) and (C 9 ) is given by 

, , C!+c 2 /t 

P V (T) = 10 


(Cll) 


When p v (l) approaches the Newtonian pressure p Q pV 2 in equation (C9), the 
contribution of the diffusion term in equation (CIO) goes to zero. When 
p v (T) exceeds pJpV 2 , the diffusion term is omitted. For the reasons dis- 
cussed following equation (C2) the computed temperature is insensitive to the 
molecular weight assigned. The value adopted is M = 0. 039^5 kg /mole . The 
most critical link in the evaluation of w(t) and the stagnation temperature 
T is the equilibrium vapor pressure law (eq. (Cll) ) , which is discussed at 
the end of this appendix. 


Allowance is made in equation (C2) for rapid evaporation or explosive 
spalling of a fraction C ms of the solid material as it reaches the surface. 
The total rate of evaporation per unit area by this means and by vaporization 
of the remaining more refractory component is given by 
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w = w + Cfo s w 8 (C12) 

Since co appears in the viscous flow equation (8), the effect of an error in 
the choice of C ms must he considered again in that connection. However , 
from the results of an investigation similar to that leading to equations (l6) 
and (17) it can he shown that an increase of C^g by a factor of 2 would lead 
to an error of less than 20 percent in the computed luminous efficiency factor. 
This finding depends partly on the choice of other parameters. With the 
parameters assigned in table I the calculations indicate that for the meteors 
analyzed the fraction of the solid evaporated is less than 25 percent of that 
removed hy viscous flow. 

This finding conflicts with the conclusion reached hy Opik (ref. 3) that 
solid stone meteoroids would evaporate rather than ahlate hy viscous flow. 

The difference arises partly from the added term l/w^ in equation (CIO) which 
often leads to a reduction in the evaporation rate helow that computed hy Opik 
based on equation (C8) alone^ The viscosity used in our calculations is also 
an order of magnitude below Opik r s value as indicated in the discussion at the 
end of this appendix. Further^ C>pik confutes the surface temperatures that 
would correspond to ablation hy vaporization alone and viscous flow alone. If 
the latter exceeds the former, he assumes that vaporization will dominate. 
However, if both ablation processes were considered simultaneously, a surface 
temperature helow either of those confuted would apply. A reduction in tem- 
perature leads to a much greater decrease in the rate of vaporization than in 
the ablation due to viscous flow. As a result, we find that efficient vapor- 
ization of stone or iron can occur only if the scale is such that ablation hy 
viscous flow is impeded hy surface tension. 

In reference 12 the skin-friction coefficient for continuum flow is deter- 
mined hy a Reynolds analogy between energy transfer and momentum transfer and 
is given hy 


Cpu = 2C Hcvu /^I (C13) 

where p 2 i is the density ratio across a normal shock and Cjj cvu is the con- 
vective heat-transfer coefficient given in equation (C 3 ). For the calcula- 
tions in this report, we have set p 21 =9* A bridging formula analogous to 
equation (C7) is 

C F = Cpue" 0 ^ + 2(1 - e " CFu / 2 ) 2 (Cl4) 

This relation yields the correct tangential force coefficient for free 
molecule flow, Cp = 2 , when Cp u becomes large. 


Viscosity and Vapor Pressure Laws 

The assignment of material properties for analysis of meteor data in 
this report is based on information contained in references 3; 15 > 21, and 22 
as well as unpublished experimental measurements. The properties of carbona- 
ceous chondrites and ordinary chondrites were considered to be of primary 
interest. Information on the properties of such material was contributed by 
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Equilibrium vapor pressure, log (Pv), N/m 2 Viscosity, log^Cfi), MKS units 
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■Tektite (ref. 15) 
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Figure 21. - Viscosity as a function of 
temper a ture. 


O Synthetic carbonaceous chrondrite 
measurements (Zimmerman) 
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ochondrite measurement (Centoianzi) 
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Figure 22. - Equilibrium vapor pressure as a 
function of temperature. 


Dean Chapman and Klaus Keil of the 
Ames Research Center. The vapor 
pressure of an enstatite achon- 
drite was measured by Frank 
Centoianzi at Ames Research Center. 
The vapor pressure of a synthetic 
carbonaceous chondrite material 
prepared by Klaus Keil was mea- 
sured by Norman Zimmerman at Ames. 
Viscosity, thermal conductivity, 
and heat capacity were measured by 
the Corning Glass Company as part 
of a larger contract monitored by 
Frank Centoianzi- 

Figure 21 shows the viscosity 
of various stones and iron as a 
function of temperature- Accord- 
ing to reference 21, log \± is 
expected to depend linearly on the 
inverse of the tenperature . The 
full line indicates the viscosity 
considered to he representative of 
carbonaceous chondrite. The small 
circles are experimental data 
points. A slope intermediate 
between that of iron and basalt 
was adopted. We have taken the 
dashed line near Opik's values for 
basalt to he representative of 
ordinary chondrites. The values 
of the parameters appearing in 
equation (Bll) based on the full 
line curve with p ms = 2,800 kg/m 3 
are C^ = 1-525X10-8 (MKS) and 
i = l4,800° K for carbonaceous 
chondrites. From the dashed line 
Cp. = 1.655X10 " 9 and i - 23,030° K 
for ordinary chondrites. The cal- 
culations in this report are based 
on the properties of carbonaceous 
chondrites (the full line in fig- 
ure 2l) . The effect on the 
results of possible errors from 
this choice is discussed in the 
Results and Discussion section. 

Figure 22 shows the equilib- 
rium vapor pressure p v as a 
function of temperature for sev- 
eral materials . From thermody- 
namic considerations log p v is 
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expected to depend linearly on the inverse of temperature with the negative 
slope proportional to the latent heat of vaporization (ref. 23). The circles 
represent unpublished experimental measurements of the vapor pressure of a syn 
thetic carbonaceous chondrite material. Because of the volatility of some com 
ponents of the mixture, it was considered advisable to allow for rapid 
evaporation of a fraction of the material and adopt a lower vapor pressure for 
the remaining more refractory part. The vapor pressure for carbonaceous chon- 
drites was therefore taken to be intermediate between the measurements and the 
values for ordinary chondrites as shown by the full line. The vapor pressure 
adopted for ordinary chondrites is shown by the dashed line. The square data 
point represents a measurement of a sample from an enstatite achondrite with 
the slope established by comparison with measurements of similar terrestrial 
materials not shown on the figure. The broken line labeled cometary meteor- 
oids is the vapor pressure used in reference 4 to compute the luminous effi- 
ciency from meteor data. The small slope was used to bring results from 
bright meteor data into agreement with the super-Schmidt data. However, the 
resulting low meteoroid mean densities for bright meteors are not now consid- 
ered realistic for such large objects on the basis of the arc-jet tests 
described in reference l4 . 

The values of the parameters c x and c 2 appearing in equation (Cll) cor- 
responding to the full line in figure 22 are c x = 10. 63 (MKS) and 
c 2 = -l6,750 o K for carbonaceous chondrites. For ordinary chondrites as rep- 
resented by the dashed line, c x = 14.215 (MKS) and c 2 = -26,700° K. The cal- 
culations in this report are based on the properties of carbonaceous 
chondrites, but the effect on the results of modifications of the equilibrium 
vapor pressure law is considered in the Results and Discussion section and 
appendix D. 
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APPENDIX D 


APPROXIMATE RELATIONS SHOWING EFFECTS OF CHANGES IN ASSIGNED PARAMETERS 


The sensitivity of the calculations to modifications of the assigned 
parameters, considered in the Results and Discussion section, has teen deter- 
mined ty approximation of equations (3)> (7)> and (8). For that purpose, it 
has teen found advantageous to utilize approximate proportionality relations 
in the place of complicated expressions* 


If the ratio T opl /T op 
trajectory for each meteor. 


in equation ( 3 ) is considered constant along the 
equation ( 3 ) can te approximated ty 


m 


„ T °pi 
T op 


m 


1 


Differentiation of equation 


(3) with respect to time yields 


(Dl) 


Equations (D2) and (Dl) lead to 


m 


T opl 

T op 


m n 



(D2) 


(D3) 


Let us consider a calculation tased on values of the assigned parameters 
different from those used to compute m. The resulting values of mass 1 % 
along the trajectory of a meteor will differ from the values of m computed 
with the parameters assigned according to table I. Equation (D3) will still 
apply, however, in the form 

ifu 

a m-L a 

The values of the ratio &i/ m i are unchanged by the reassignment of param- 
eters since they arise directly from the input quantities (h, V, V, m.) taken 
from reference 5* Therefore, dividing the last relation by equation (D3) 
yields 

m m 

This can be expressed as the proportionality relation 

m ex m (D4) 


which applies with respect to variations of the assigned parameters. That is, 
the computed value of m on a point of the trajectory of a particular meteor 
will change by the same fractional amount as the computed value of mass 
changes due to a reassignment of parameters. 
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For small variations of p ms or pjuf from their assigned values, the 
factor 

Pm ~ Pmf 
y ~ Pms " Pmf 

appearing in equation (7) is approximately proportional to a product of powers 
of Pm> Pms' and Pmf 


Pm 

Pms 


p mf 

p mf 


^PmPmsPmf 


where C is a constant. Variations of % from its computed value must be 
considered because of the effect of changes in p^g or Pm-p in other equations. 
The operation 8ln(y)/8p m performed on both sides of the last equation and 
rearrangement leads to 

_ Pm 
Pm “ Pmf 

Evaluating the right side at the typical vsflues % = 0.4 g/em 3 and 

= 0.1 g/cm 3 yields a = 4/3. Similarly, partial logarithmic differentia- 
tion with respect to Pujg and pjjjf with p ms = 2.8 g/cm 3 leads to b = -28/27 
and c = -8/27. Since the initial confuted value of Pm can differ from 0.4, 
the exponents b and c have been simplified to -1 and -1/3. Therefore, the 
following approximate proportionality is arrived at 


Pm ~ Pmf 
Pms - Pmf 


4/ 3 -1-1/3 


Tn 


ms 


p mf 


(D5) 


For changes in p m , p ms , or p^ by a factor of 2 from their initial values, 
the error in this approximation is less than 20 percent, although each of the 
quantities on the two sides changes by as much as a factor of 3, Gross 
inaccuracies arise only when approaches such that the left side 

changes by an order of magnitude. Equation (D5) applies with respect to varia- 
tion of any assigned parameter. The factors p^g and P^f /3 could be dropped 
out of the right side of equation (D5) if we wished to consider the effects of 
variations of other assigned parameters while holding p ms and p^ fixed. In 
contrast, the factor p^ /3 could not be dropped for any case because p m is a 
computed quantity rather than an assigned quantity. 


From the machine calculations it has been ascertained that the term in 
equation (7) containing is ordinarily small compared to the first term 

containing m. In that case the approximate proportionality relation 

-2 -1 1/9 2/3 -1/9 1/3 , 

oj s « A rs S F p m p ms m U 

can be derived from equation (7) with the aid of equations (d 4) and (D5) . 

Equation (8) leads to the proportionality 


TH 1 


S F 


1/2 1/3 -1/3,2 


m 


' “)C; 1 (PV 2 ) 


(D7) 
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where 


F = 


ikf 

%s) 


Pms 

^(T) 


(T - T f ) 2 G(T) 


(D 8 ) 


The factor (pV 2 ) -1 need not he retained for variations of the type under con- 
sideration since it depends only on the input data (h, V, V, m^,) and would not 
he affected by changes in any assigned parameter* Also, to is ordinarily 
small compared to w s so that the relation 




i 0 -i/2 -1/3 1/3. 




m 


p m fc f 


can he derived from equation (D7). 


yields 


Substituting equation (D 6 ) into the last equation and rearranging it 
.s 


m 


i 4/3 o2„*xJ /s « aZj&Lc-fsI 72 


Further rearrangement leads to the relation 


.- 3/4 9/2 15/8 - 3 / 2 . 1/4 3/4 3/4 
m oc A rn Apg Sp p ms p^ Cp F 


(D9) 


All quantities appearing in this expression are assigned parameters 
except m, Cp, and F which are computed, hut also depend on additional 
assigned parameters. However, it can he shown that Cp and F are rarely sen- 
sitive to changes in the assigned parameters contained explicitly in equa- 
tion (D9) . Therefore, Cp and F can he dropped from equation (D9) for pur- 
poses of considering such variations. The effect of changes in the assigned 
parameters contained in Cp and F can he considered separately hy omitting 
all factors except Cp and F. In other words, equation (D9) can he factored 
into the two relations 


- 3/4 9/2 15/8 - 3/2 1/4 

m cc A rn A rs op Pms Pmf 


(DIO) 


and 


3/ 43/4 
m ^ Cp F 


(Dll) 


hy virtue of an approximate independence of the two factors with respect to 
the types of variations of assigned parameters under consideration. 


The nature of the approximation can he understood hy considering a varia- 
tion of A rn . According to equation (DIO), this would change the computed 
values of m, which would change the radii r n and r in equation (C 3 ), so 
that Chcvu anc ^ Cp evaluated hy equations (C13) and (Cl4) would he altered. 
Therefore, strictly, Cp depends on m and should he retained in equa- 
tion (D9) which would invalidate equation (DIO) . However, equation (C 3 ) shows 
that Chcvu depends on m 1/6 (since r n m 1/3 ) if we exclude the few cases 
where the exponential term causes a stronger coupling. Therefore, an accurate 
evaluation with Cp included in equation (DIO) would differ only slightly 
from the result obtained hy excluding j it , except possibly for a few points on 
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the trajectories of a small number of meteors. The statistical procedure 
would tend to suppress the effect of the exceptional cases in the computed 
averages. 

The quantity F is defined in equation (d 8). For a particular meteor, 

F depends only on the stagnation-point temperature T and the assigned 
parameters contained in the material property laws. In other words, F 
depends on the other computed quantities such as m only through their effect 
on T. But it can he seen from equation (C2) and the discussion thereafter 
that this dependence is weak. Therefore, omitting F from equation (DIO) is 
not expected to cause serious error in the consideration of the effects of 
variations of the assigned parameters A rn , A rs , Sp, p^, and Pmf • IT we wish 
to consider the effect of variation of only one of the assigned parameters in 
equation (DIO), the others can be omitted from that equation since they are to 
be held fixed in the variation. For that reason, equation (Dll) follows from 
equation (D9) under the condition that the assigned parameters A r nj A rs , Sp, 
p ms , and are to be held fixed. 

The derivation of equations (DIO) and (Dll) is based on the assumption 
that the stagnation tenperature T is approximately independent of variations 
in the assigned values of geometric parameters. This means that T is 
assumed to depend only on input quantities and assigned parameters appearing 
in the aerodynamic and material property laws. That dependence will be con- 
sidered explicitly so that equation (Dll) can be interpreted in terms of varia- 
tions of the assigned parameters. However, it is worth noting at this point 
that, by virtue of the approximations employed, the values of mass computed 
have been found to depend only on the assigned parameters contained in the 
ablation model. Within the accuracy of the approximations of this appendix, 
the values of mass do not depend on the luminosity and dynamic equations. 

The relation 

m cc -r~p (D12) 

follows from equation (Dl). Substituting this into equation (DIO) leads to 
equation (l6) of the Results and Discussion section. Again, it is worth 
noting that within the accuracy of the approximations of this appendix neither 
the mass nor the luminous efficiency factor depends on the dynamic equation. 
This means that values of the luminous efficiency factor and mass computed by 
the machine program are expected to be insensitive to errors in the accelera- 
tions that may exist in the reduced meteor data used as input. The effect of 
such errors, if present, would be confined largely to the computed values of 
mean meteoroid density. 

In order to determine the effect on the mass and the luminous efficiency 
factor of changes of the assigned parameters in the aerodynamic and material 
property laws, a further analysis of the quantity F defined in equation (d 8) 
is needed. The first step in this procedure is to determine how changing 
these parameters affects the computed stagnation-point temperature. 

The term eoT 4 in equation (C2) is ordinarily small compared to the 
other terms. If it is omitted, the temperature can be evaluated in two steps: 
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(l) The vapor pressure p vc is determined by equations (C2) to (CIO) • The 
subscript c is intended to indicate that the vapor pressure under considera- 
tion is computed from the aerodynamic relations and does not depend on the 
assignment of parameters in the equilibrium vapor pressure law. ( 2 ) The value 
of p vc from the first step is used to determine the stagnation-point tempera- 
ture T from equation (Cll). For that purpose, equation (Cll) can be 
rearranged in the form 


C 2 

C 1 rp “ PvC 


Let us consider the effect on T of variations in c x or p vc and let 
the subscript b refer to the base values corresponding to the assigned 
values in table I. The relation 


' 2 b 

T 


= c 


lb 


Ti 


2b - Pvc 
108 Pvt 


is obtained by combining two equations of the above form and indicates the 
effect on T of changes in c ± or pvc • Rearrangement leads to the equation 

( T b ~ T) ^ Pvc 

c 2b TT ~ c lb c i^" ^ 

1 b 1 p vb 

Since T will not differ greatly from T^, a typical value of 2,400° K can be 
used for both T^ and T in the denominator on the left. Upon insertion of 
the value c^ = - 16,750 from table I as well, the relation 


can be arrived at. 


T - T. = 3^0 


xb 


c-l + log 


PycN 

Pvb/ 


The last expression indicates that if the equilibrium vapor pressure law 
(eq. (Cll)) is modified so that the level of vapor pressure at a given tem- 
perature is higher by a factor of 10 (c ± = c^ + l) , the computed temperatures 
will all be decreased 3^0° K. The expression also indicates that if the value 
of p vc computed from the aerodynamic relations is increased by a factor of 
10 the computed temperature will increase 3^0° K. It can be seen that reason- 
able reassignments of the parameters in the aerodynamic laws, equations (C2) 
to (CIO), lead to changes in p vc by a factor of order 2 rather than 10. 

This would change the computed temperatures about 100° K. Changes in p vc 
due to alterations of other computed quantities such as m or r are 
ordinarily much less than a factor of 2. 

It will be seen that a change of 100° K leads to a modification of about 
30 percent in the computed values of the luminous efficiency factor. There- 
fore, we shall drop the consideration of effects of modifications other than 
the level of the equilibrium vapor pressure law (the value of cj and rewrite 
the last expression as 

P u 

T - T b = 3^0 log — (D13) 

Pv 
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where - c has been replaced by logCp^/p-y-) for notational convenience. 

This relation will he used to evaluate the temperature appearing in the 
quantity F. 

Substituting equation (Bll) into (D 8 ) yields 


F = 


(^) 2 C^ 1 e~ 1//T (T - T f ) 2 G(T) 


The function G(T) defined in equation (B13) can he approximated hy 


so that 




C e ~ l/T(T ' 


) 2 ' 


©‘ 



The dependence of F 
he approximated hy 


on all quantities except i 

T-Tb 



in the above expression can 


for small deviations of the assigned parameters or T from their hase values. 
The quantities h and T^ are constants to he evaluated and C is a propor- 
tionality constant. If partial derivatives with respect to Tp and T of the 
logarithms of the last two expressions for F are equated and evaluated at 
the typical values T = 2 ,400° K, Tp = 1,800° K, i = l4,800° K, one finds 
h = 0.25 and T ^ = 292° K. Substituting these into the last expression for F 
yields the proportionality relation 


F CC (^ C -V* e 292 ( Dl4 ) 

Using this expression in equations (Dll) and (D 12 ) shows that a change in tem- 
perature T - T^ = 100° K leads to less than a 30-percent change in T Q p - 

Substituting equation (DI 3 ) into (Dl4) replaces the temperature dependent 
exponential factor with the factor (pyb/Pv) ° # 505 • Since the typical value of 
T used in the derivation is not exact, the exponent O. 5 O 5 can he replaced by 
l/2. Also, the hase vapor pressure p^ corresponding to the assigned param- 
eters in table I can he dropped from the proportionality relation since it 
does not depend on any reassigned parameters. The resulting proportionality 
relation for F is 

58 


i 



F 


'k \ 2 -11/4 -i/; 
— 1 C M T- -n 


£ f J t P 


Substituting equation (D15) into (Dll) yields 


m oc Cp 


/ 4/j L \ 

\£f/ 


3/2 




3/4 3/16 -3/8 


t: 


■V 


(D 15 ) 


(di 6 ) 


This expression contains only assigned parameters except for the tangential 
force coefficient Cp and p v - As discussed earlier, we shall neglect changes 
in Cp due to changes in computed quantities; Cp is included to indicate the 
effect of possible errors in the assigned parameters that affect its value 
such as the constant C in equation ( 03 ) * The symbol p v has been used here 
in place of 10 Cl , which appears in the equilibrium vapor pressure law 
(eq. (Oil)) and is hence an assigned parameter. Substituting equation (D12) 
into (Dl6) leads to equation (lT) of the text. 


Additional approximate analyses and machine calculations have been made 
to determine the effects of changes in the parameters c 2 and i appearing in 
the viscosity and vapor pressure laws. Such variations were found to be 
approximately equivalent to changes in the parameters Cjj. and p v considered 
above except for an additional constant shift in the stagnation-point tenpera- 
ture . Machine calculations were also made with different values of the 
parameters appearing in equations (C2) through (C9) * such as e, £ v , and C ms . 
For changes by as much as a factor of 2, the effects on m and T 0 p were less 
than 30 percent as anticipated in the discussion preceding equation (D13)* 


Results from the approximations of this appendix are conpared with 
results from machine calculations in the Results and Discussion section. The 
comparisons show that the approximations are reasonably accurate in the pre- 
diction of the effects of reassignment of all parameters except the froth den- 
sity Pxnf* * Considerable insight into the problem under investigation is 
thereby afforded by the approximations as summarized in equations (l6) and 
(l7) of the text. How reassigning parameters affects computed quantities, 
such as mass, mean meteoroid density, or stagnation-point surface temperature, 
can be deduced from formulas contained in this appendix. 
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APPENDIX E 


INTENSITY LAG EFFECT 


Relationship Between T 0 p C and T 0 p 
Equation ( 19 ) can he rewritten in terms of the known quantities 

noo j 

“i(t) = « (El) 

T Opl^t v 

where T 0 pj is the assumed constant value of luminous efficiency factor upon 
which the values of mass mi published in reference 5 are based. With the 
assumption that T 0 p C is constant along the trajectory and upon substitution 
of equation (El), equation ( 19 ) can be written 


m(t) = ° Pl {m^t) + Btm^t + 0) - m 1 (t)]} 
T opc 


For small Q , m^t + 0 ) can he approximated by 

m-j/t + 0) = m-j^Ct) + 0m 1 (t) 
and substituting this into equation (E2) yields 


, , T opl 
m(t) = m 


T opc 


^(l + B0 


(E2) 


(E3) 


(E4) 


(E5) 


Differentiating and combining equations (l8) and (El) leads to the 
relation 

Top 

= m 

1 T opl 

If T 0 p in equation (l8) is taken to be approximately constant along a tra- 
jectory so that T 0 p can "b e taken outside the integral, combining the result 
with equation (El) yields 


mi = 


Top 


m 


T opl 

Substituting equations (E 5 ) and (E6) into (e 4) leads to the approximate 
relation 


(E6) 


opc ‘op 

which is equation (20) of the text. 


= r_( 1 + BO - 
m 
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lag Time 


An equation for the rate of ablation of froth particles in the wake can 
be written 


2_ duif 2_ Cj^pV 3 Ap 

mf dt 2 £ m-f 


(ET) 


Substituting equations (4) and ( 5 ) leads to 


where 


in 


-2/ 3 


dmf 

dt" 


= -K 


K f = 


(l/2)C H pV 3 S F 



The integration of equation (e 8) with Kf 


considered constant yields 


(E8) 

(E9) 


1/3 1/3 1 v /, , \ 

“fo “ “f = 3 K f( t " t o) 


(E10) 


This relation indicates that the froth particle would ablate to nothing in a 
time 


Substituting equation (E9) yields 


C^pV 3 Sp 

which is equation (2l) of the text. 

According to reference 3; chapter 6, under air pressure , liquefied 
meteoric material will break up into drops of radius 




3m. 


0 = 


1/3 

To 


% 


2 /3 l/3 






(Ell) 


r 


l 



(E12) 


where t s is the surface tension. In reference 15 > the numerical value 2.6 
is used in place of 4.8 in the above expression to indicate incipient disrup- 
tion of a tektite drop . A value of the surface tension 


t g = O .36 kg/sec 2 


is used for stone or pumice in reference 24. Ablation can occur only by 
vaporization for drops of the above size or smaller. Therefore, the specific 
heat of ablation £ in equation (2l) should be taken equal to £ v for such 
drops . 


61 


Substituting 


4 3 

m l = 3 1 %f r ? 

and the shape factor for a sphere 



into equation (2l) leads to 


Evaluating this with £ v = 8 . 5 IXIO 6 , = 100 kg/m 3 , t s = O .36 kg/sec 2 ,, 

V = 2X10 4 m/sec, Cjj = 1 yields 

9 Z = 0.0235x(l0 8 /pV 3 ) 2 

which is equation ( 25 ) of the text. 

At an early part of the trajectory of meteor 4464 of reference 5 the 
impact pressure and velocity were 

pV 2 = 3^00 N/m 2 

Y = 2X10 4 m/sec 

Under these conditions and with the above value of surface tension, 
equation (E12) indicates 

r^ = 5*1X10~ 4 meter 

The corresponding values of mass m^ and Gi are indicated in the text. 

For solidified froth particles, an estimate of the initial mass m.f 0 is 
needed to determine Of from equation (2l) . The faintest meteors described 
in reference 20 show fluctuations in the light curves at a rate of about 15 
fluctuations per second. According to the method of this report, meteor 4464 
loses mass at about the rate of 4.5X10~ 4 kg/sec. If all the ablation is 
attributed to sloughing of froth particles, the mass per particle would be 
about 

mf Q = 3X10" 5 kg 

Evaluation of equation (2l) with = 3X10 5 kg and £ = £ v yields 

0 = 0.42 sec (vaporization dominant) 

which is somewhat high compared to the value of the product B0 = 0.2 second 
obtained for meteor 4464. Use of the viscous flow equations derived in appen- 
dix B leads to the conclusion that, for particles of the size and composition 
under consideration, ablation will occur largely by liquid run off rather than 
by evaporation. In that case the specific heat of ablation £ in 
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equation (2l) would be closer in value to £f than The resulting 

computed lag time is 

Of = 0.092 sec (liquid ablation dominant) 

Fragmentation of the froth particles may also play a role in producing a 
higher rate of vaporization of the particles in the wake than that occurring 
directly from the parent meteoroid. 

The dependence of Of on pV 3 is found by substituting £ = 1.884x10 s , 
Pmf ” 100 kg/m 3 , m fo = 3X10" 5 kg, = 1, S F = (9Jt/l6) 1/3 in equation (2l) to 
obtain 

Of = 0.063X10 8 /pV 3 
which is equation (24) of the text. 


Lag Distance 


A dynamic equation for the froth particles similar to equation ( 2 ) can be 
written 

- - i ^ ^ 

Since the combination CpA/m^ is larger for the froth particles than for the 
parent meteoroid, they would be expected to lag behind once they were free of 
the protection of the stream in the immediate wake. According to W. A. Page 
of Ames Research Center, thin plastic disks have been observed to overtake and 
bump models of smaller CpA/m in free -flight tests. Thus it is possible that 
during deceleration an accumulation of froth builds up behind the parent 
meteoroid. However, we assume that froth particles can also escape into unpro- 
tected regions of the wake due to bumping or other dispersive influences. In 
that case, they will accelerate away from the parent meteoroid approximately 
at the rate 

Ai f - \ C D pT 2 2 (El*) 


If AVp is taken to be constant, the separation distance is 

AL = ~ AV^»(t - t Q ) 

and at a time t - t Q = 0 will amount to 

AL = | AV f e 2 

Substituting the above expression for AVf and using equation (E9) to 
remove pV 2 leads to 

1 c d £q 
^ = 2 1/3 c H v 

m fo H 
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Substituting equation (Ell) to remove 


Kf yields 


AL 


3 

- 2 C H V 


•which is equation (22) of the text. An expression for AL with a numerical 
factor of 3 rather than 3/2 can be derived by considering the variability of 
Af/mf in equation (El4) rather than setting it equal to Sp/(Pmf 3 ^fo 3 ) • How- 
ever , this value would be reduced if the smaller values of pV 2 at the early 
part of the ablation process were taken into account. Therefore, we have used 
the above numerical value of 3/2 for estimating the lag distance. 
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TABLE I.- ASSUMED VALUES OE CONSTANTS AND PARAMETERS 3 - 


Quantity 


Value , 
MKS 


Name 


A rn 

1.0 

Shape parameter 

A rs 

0.874 

Shape parameter 

C 

0 . 943 X 10 -4 

Heat -transfer parameter 

Ci 

10.63 

Parameter in vapor pressure law 

C 2 

- 16,750 

Parameter in vapor pressure law 

Cms 

0.15 

Parameter in energy equation 

Cp 

1.525x10" 8 

Parameter in viscosity law 

i 

l 4 , 8 oo 

Parameter in viscosity law 

k 

2.0 

Thermal conductivity 

M 

0.03985 

Molecular weight 

Pv 

0.15 

Heat-transfer parameter 

Sp 

1.919 

Shape factor 

T f 

1785 

Fusion temperature 

a 

0.3 

Heat-transfer parameter 

e 

0.6 

Radiative emissivity 

P 21 

9 

Density ratio for normal shock 

Pmf 

100 

Froth density 

Pms 

2,800 

Solid density 

a i 

0.1 

Heat-transfer parameter 

£f 

1 . 884 xiO s 

Specific heat of fusion 

£y 

8.51x10 s 

Specific heat of vaporization 


a Material properties correspond to carbonaceous chondrite • 


TABLE II.- TRAIL NUMBERS IN 40-METEOR SAMPLE 


3,076 a 

8,447 

8,945 

9,880 

10,240 

11,825 

12,399 


4 , 464 a 
8,469 

8,990 

9,888 a 

10,358 

11,856 

12,504 a 


6 , 9 4 9 a 

8,726 

9 , 015 a 

9,900 

10,384 a 

11,973 

12,577 


7,l6l a 
8,766 a 
9,030 
9,917 
10 ,4l4 a 
12,342 
12,714 


7,2l6 a 

8,88i a 

9,170 

10,070 a 

10,439 

12,381 


7,272 a 

8,891 

9,4l6 a 

10,094 

10,447 

12,363 


Q - ■ - I 

Meteors for which intensity lag times and corrected values 
of luminous efficiency factor were computed. 


A- 3004 


NASA-Langley, 1968 30 
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